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Abstract 

The extragalactic background light (EBL) is one of the fundamental ob- 
servational quantities in cosmology. All energy releases from resolved and 
unresolved extragalactic sources, and the light from any truly diffuse back- 
ground, excluding the cosmic microwave background (CMB), contribute to 
its intensity and spectral energy distribution. It therefore plays a crucial 
role in cosmological tests for the formation and evolution of stellar objects 
and galaxies, and for setting limits on exotic energy releases in the universe. 
The EBL also plays an important role in the propagation of very high en- 
ergy 7— rays which are attenuated en route to Earth by pair producing 7 — 7 
interactions with the EBL and CMB. The EBL affects the spectrum of the 
sources, predominantly blazars, in the ~ 10 GeV to 10 TeV energy regime. 
Knowledge of the EBL intensity and spectrum will allow the determination 
of the intrinsic blazar spectrum in a crucial energy regime that can be used 
to test particle acceleration mechanisms and VHE 7— ray production models. 
Conversely, knowledge of the intrinsic 7— ray spectrum and the detection of 
blazars at increasingly higher redshifts will set strong limits on the EBL and 
its evolution. This paper reviews the latest developments in the determina- 
tion of the EBL and its impact on the current understanding of the origin 
and production mechanisms of 7— rays in blazars, and on energy releases in 
the universe. The review concludes with a summary and future directions 
in Cherenkov Telescope Array techniques and in infrared ground-based and 
space observatories that will greatly improve our knowledge of the EBL and 
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the origin and production of very high energy 7— rays. 

Keywords: extragalactic background light, cosmic infrared background, 
cosmology, dark matter, galaxy evolution, gamma-ray astronomy, 

GeV /TeV sources, blazars, gamma-ray opacity 


1. Introduction 

The extragalactic background light (EBL), dehned here as the emission in 
the 0.1 to 1000 /rm wavelength region, is one of the fundamental observational 
quantities in cosmology. It comprises the integrated light from resolved and 
unresolved extragalactic sources, and the light from any truly diffuse back- 
ground, excluding the cosmic microwave background (CMB). It is therefore 
the repository of all energy released by nuclear and gravitational processes 
since the epoch of recombination. A significant fraction of this radiation is 
shifted by cosmic expansion and by absorption and reradiation by dust into 
infrared (IR) wavelengths. Consequently, its intensity and spectral shape 
hold key information about the formation and evolution of galaxies and their 
stellar and interstellar contents throughout cosmic history. A strict lower 
limit on the EBL intensity is provided by the integrated light from resolved 
galaxies, hereafter referred to as the integrated galaxy light (IGL). 

The EBL plays also an important role in the propagation of high energy 
7— ray rays that are predominantly emitted by blazars, a subgroup of active 
galaxies hosting active galactic nuclei (AGN), whose relativistic jet is pointed 
towards the Earth. High energy photons emitted by blazars are attenuated 
by photon-photon interactions with the EBL, a process that can be used to 
set important limits on both, the intrinsic spectra of blazars and the intensity 
of the EBL in select energy and wavelength regions where these interactions 
are most prominent. 

The EBL is intimately connected to the diffuse X-ray, radio, and super- 
nova neutrino backgrounds. Deep X-ray surveys have resolved the X-ray 
background into point sources, most of which are dust enshrouded AGNs 
(Mushotzky et ah, 2000). Up to 90% of the X-ray energy produced in indi- 
vidual AGN can be degraded and reradiated predominantly at mid-IR wave- 
lengths (e.g. Franceschini et ah, 2002; Ballantyne et ah, 2006). Gonsequently, 
the X-ray background can be used to predict the EBL intensity at at these 
wavelengths. Gurrent estimates show that about 15% of the 24 /rm EBL in- 
tensity is powered by AGN activity (Treister et ah, 2006; Soifer et ah, 2008, 
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and references therein). Conversely, the connection between mid-IR bright 
sources and AGN can be used to estimate the contribution of obscured AGN 
to the X-ray background (Gandhi & Fabian, 2003; Soifer et ah, 2008, and 
references therein). 

Massive stars that power the IR emission also emit radio free- free emission 
during the main sequence phase, and radio synchrotron emission during the 
supernova remnant phase of their evolution. The IR emission from star- 
forming galaxies is therefore correlated with the radio emission (Lisenfeld 
et ah, 1996; Condon et ah, 1991). This correlation can be used to estimate 
the contribution of star-forming galaxies to the cosmic radio background 
(Haarsma & Partridge, 1998; Dwek & Barker, 2002; Ponente et ah, 2011). 

Most of the EBL intensity is powered by massive stars that end their 
life as core collapse supernovae. The total EBL intensity can therefore be 
used to derive an estimate of the supernova rate and the resulting flux of 
supernova neutrinos (Horiuchi et ah, 2009; Beacom, 2010). The detectability 
of these neutrinos can be greatly enhanced by the proposed introduction 
of gadolinium in existing large water Cherenkov detectors (such as Super- 
Kamiokande) (Beacom & Vagins, 2004). Gadolinium has a very high capture 
cross section for neutrons generated in Pe + P — t reactions, and can be 

introduced in the form of soluble trichloride (GdCls). Following the neutron 
capture, the Gd emits an 8 MeV 7— ray which produces relativistic electrons 
by Compton scattering. The Cherenkov radiation from these electrons is 
more easily detected than that produced in the cascade of the 2.2 MeV 
7— ray generated by the capture of neutrons by free protons. 

Several reviews have appeared in the literature, presenting a historical 
overview of the importance of the EBL, early estimates of its intensity, the 
quests for its detection, and its many astrophysical implications (Hauser & 
Dwek, 2001; Kashlinsky, 2005; Lagache et ah, 2005). Since these reviews 
were written signihcant advances have been made in studies of the EBL with 
the launch of UV {Galex) and IR space observatories {Spitzer, Herschel, and 
Akari). These observatories, together with ground-based telescopes, such as 
2MASS, have provided new limits on the EBL ranging from UV to submil- 
limeter wavelengths. Deeper galaxy number counts and new data analysis 
techniques of stacking astronomical images have narrowed the gap between 
the contribution of resolved galaxies and the true intensity of the EBL. 

The Fermi Gamma-ray Space Telescope, operating between 200 MeV and 
300 GeV, and ground-based air Cherenkov detectors (H.E.S.S., MAGIC, and 
VERITAS) operating in the ~ 50 GeV to 100 TeV range have broadened 
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the energy window for the studies of 7— ray sources. These advances have 
led to the detection of new GeV and TeV 7— ray sources and provided new 
data for determining their intrinsic spectra. Reviews of these subjects were 
presented by Weekes (2008) and Hinton & Hofmann (2009). More recently, 
Dermer (2012) presented a review of the Fermi catalog of 7— ray sources and 
the physics of the production of relativistic particles and 7— rays from these 
sources. Table 1 presents a glossary to the acronyms of the observatories and 
instruments referred to in this review. 

These developments provide the main impetus for this review. We hrst 
present, in §2, the basic formulae describing the attenuation of photons by 
pair producing interactions with other photons. We then show how this at- 
tenuation will affect 7 rays traversing a radiation held characterized hrst by 
a pure black body, representing the stellar emission component of the EBL, 
and then by a more realistic EBL that includes the dust emission component. 
This attenuation can, in principle, be used to determine the intensity of the 
attenuating radiation held if the intrinsic source spectrum is known. In §3 we 
survey the type of 7— ray sources that are used in these studies, their spec- 
tral characteristics, the physical mechanisms for generating their spectra, and 
constraints on their spectral shape imposed by general physical principles. In 
§4 we summarize measurements and limits on the EBL intensity determined 
by direct measurements and by adding the light from resolved galaxies. Mod- 
els for the EBL intensity and its evolution with redshift are summarized in 
§5. In §6 we summarize the constraints on the EBL intensity derived from 
7— ray observations of blazars, emphasizing the diherent assumptions made 
on the intrinsic blazar spectra to derive these limits. EBL models predict 
the 7— ray opacity of the universe at different energies, and in §7 we compare 
these model predictions with blazar observation. Throughout this review it 
was tacitly assumed that the production of 7— rays takes place exclusively in 
the sources. In §8 we consider alternative scenarios of 7— ray production that 
could have important implications for EBL limits, namely, that a signihcant 
fraction of the observed 7— rays could be produced en route to Earth. The 
role of the EBL in setting limits on exotic energy releases in the universe 
in briefly discussed in §9. A summary and future prospects for the fields of 
7— ray and EBL research is given in §10. 
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2. The EBL and the Attenuation of Gamma-Ray Photons 

2.1. The EBL 

The differential speciffc ffux at wavelength Aq, dFi,(Ao), received from 
radiative sources within a comoving volume element dV'c(^) at redshift at 
wavelength A is given by (e.g. Mo et ah, 2010); 


dF,(Ao) 


Cu{\,z) dVcjz) 
47T d.L{zy 


( 1 ) 


where z) is the comoving speciffc luminosity density of the sources, d,L 
is their luminosity distance, and the (1+z) factor arises from the decrease in 
energy of the emitted photons due to the redshift, and Aq = (1 + z)\. 

The speciffc comoving intensity of the EBL per unit solid angle, 5f2, at 
redshift Zq and wavelength Aq is given by an integral over all energy releases 
over cosmic history: 


d;/(Ao, Zq) 



Cu{X,z) dVc{z) 
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dz 


where c\dt/dz\ is given by (e.g. Mo et ah, 2010): 

f? ^ JL 
{l + z)E{zy ^ Ho 

+ fim(i + + f2fc(i + ^y + iia] 

= [(1 + zy (LlmZ + 1) — z{2 + z) IIa] 

= + zy + Ha] ^ 


dt 

dz 


( 2 ) 


( 3 ) 

( 4 ) 


Ho is the Hubble constant, and and Ha are the dimensionless 

density parameters of the radiation, matter, the curvature, and the cosmo- 
logical constant A, obeying the relation: H/? -|- H^ + H^ -|- Ha = 1. The second 
expression for E{z) is for a matter dominated {Qr << 1) universe, and the 
third is for one that is matter dominated and ffat (H^ = 0). In the concor- 
dance cosmology model: Ho = 70 km s“^ Mpc“^; H^ = 0.27, and Ha = 0.73 
(Hinshaw et ah, 2009). 
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2.2. Gamma-ray attenuation by pair production 

The interaction between two photons with energies E.y and will lead 
to the creation of a particle anti-particle pair when the total 7— ray energy 
in the center of momentum of the system exceeds the rest frame energy of 
the two particles. The threshold for the creation of an e^-l-e“ pair is given 
by: 


P,, Z'^ 


2 (rUe 
^7 (1 - /^) 


(5) 



Figure 1: Schematic illustration of the 7 — 7 pair production reaction, showing the defini- 
tion of the angle 6 between the interacting photons. 


where p = cos6', and 9 is the angle between the two photons, as illustrated 
in Figure 1. 

The cross-section for the 7 — 7 interaction is given by: 


3(X 

e, p, z) = (1 - (3'^) 


2/3 {(3^ - 2) + (3 - (3^) In 


1 + / 3 ) 
(1-/9) 


( 6 ) 


where 



(7) 


Figure 2 (left panel) depicts the cross section as a function of (3. The 
cross section peaks at a value oi (3 = 0.70, providing a relation between the 
energies E..^ and e (or wavelength A) at the peak, given by: 


E^iTeV) 


1.07 

e{eV) (1 -p) 


0.86 A(pm) 

(1 -p) 


( 8 ) 


The right panel of the hgure depicts the cross section as a function of 6 = 
2{mc^Y / E.yC for different values of the angle 6. When the photons are moving 
in the same direction (0 = 0), the cross section collapses to a delta-function 
at 6 = 0, and the energy threshold becomes inhnite. 
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Figure 2: The cross section for the 7 — 7 interaction. Left panel: its dependence on /3 
[eq. ( 7 )]; Right panel: its dependence on b for different angles of incidence. 

2.3. The Attenuation of 'j— rays from Cosmological Sources 

En route to Earth, 7— rays from cosmological sources have to pass through 
the radiation held of the EBL, resulting in their attenuation by pair producing 
interactions. The optical depth of a 7— ray photon at an observed energy 
emitted by a source at redshift z due to this process is given by: 

r-,j{E.y,z) = dz' ^ J demC, z'){l + z'f , z') ( 9 ) 

where n^{e, z) = dn{e, z) / de is the specihc comoving number density (cm“^ eV“^) 
of background photons with energy e at redshift z, and the (1 + z)^ term 
represents its conversion to a proper number density. The pair-production 
threshold energy is = 2{meC^)'^ / E^{1 — /i)(l -|- z), where the (1 -f z) fac- 
tor takes into account that the observed 7— ray photon had a higher energy 
at the redshift of the interaction. The parameter /3' = (1 — and 

di/dz = c\dt/dz\, where i is the proper distance. 

Calculating the EBL opacity to 7— rays from cosmological distant sources 
requires knowledge of the evolution of the comoving specihc photon number 
density n^{e, z) as a function of redshift. The specihc number density of 
photons with energy e at redshift 7 is related to the specihc EBL intensity 
at a given redshift 7 by: 

47T 

e^ne{e,z) = — u C{u, z) (10) 

c 

= 2.62 X 10"V4(z/,7) 

where e = hu, 7) is given by eq. (2), and the coefficient in the second 
line was calculated for e in eV, in cm“^ eV“^, and u C in nW m~^ sr“^. 
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Finally, we point out that the 7 — 7 cross section is wide, so that in 
calculating the 7— ray opacity, strong variations in the EBL spectrum are 
smoothed out over a wide range of 7— ray energies. The EBL intensity at a 
given wavelength is therefore effecting over a wide range of 7— ray energies 
around the peak given by eq. (8). 

2.4. A Simple Example: An EBL given by a diluted blackbody spectrum 

Of particular interest is the behavior of for a background radiation 
field that is represented by a diluted blackbody. Eigure 3 (upper left panel) 
depicts a local EBL characterized by a Planck function, normalized to an 
intensity of 10 nW m~^ sr“^ at 1 pm. The upper right panel of the figure 
depicts the photon number density. The bottom left panel shows the 7— ray 
opacity at redshift = 0.2, assuming a non-evolving EBL, and the right 
panel shows the source attenuation as a function of 7— ray energies. Also 
shown in the hgure are the energy regimes in which substantial changes in 
the slope of the opacity occur (dashed lines). 

The rapid rise in the EBL spectrum between 0.5 and 1 pm results in a rise 
of the 7— ray opacity, and the onset of substantial source attenuation in the 
10 to 500 GeV energy region. This sudden increase in the GeV attenuation 
creates a break, TceVi in the spectrum, dehned as the difference in power law 
index between the unattenuated and the attenuated region of the spectrum 
(see Figure 5 in this paper). At higher 7— ray energies, the spectrum of a 
blazar characterized by an intrinsic power law will exhibit a second spectral 
break around ~ 1 TeV. For an evolving EBL, the magnitude and location 
of this spectral break are expected to evolve with redshift. The substantial 
decrease in the attenuation at a few TeV is a consequence of the particular 
choice of the EBL spectrum , which decreases rapidly at wavelengths beyond 
~ 2 pm. 

2.5. A More Realistic Example: An EBL that includes dust emission 

Figure 4 depicts a more realistic presentation of the current EBL spectrum 

(left panel) and the 7— ray opacity for different redshifts (right panel), taken 
from model calculations of Finke et al. (2010). At wavelengths short wards 
of ~ 5 pm the spectrum represents the stellar and AGN contributions to the 
EBL. At longer wavelengths the spectrum represents the AGN and starlight 
energy that was absorbed and reradiated by the dust. The right panel shows 
the energy dependence of the 7— ray opacity for sources at different redshifts. 
The opacity calculations took into account the evolution of the EBL with 



z = 0 






E^(TeV) E/TeV) 


Figure 3: Top left; A diluted black body representation of the stellar emission component 
of the EBL; Top right: The corresponding proper photon number density versus energy; 
Bottom left: The 7— ray opacity versus energy, Ej-, Bottom right: The 7— ray attenu- 
ation. The figure illustrates the dramatic change in the attenuation at the 7— ray energy 
that corresponds to the wavelength at which the slope of the EBL spectrum changes. The 
different slopes are depicted as dashed lines in the figure. 


redshift. The figure illustrates the relation between the EBL spectrum and 
the energy dependence of the 7— ray opacity. The initial rise of the EBL 
intensity at UV-optical wavelengths causes an increase in the 7 — 7 opacity 
between 10 and 500 GeV. The decline in the EBL intensity between ~ 1 and 
15 /im causes to rise less rapidly between 1 and 10 TeV. The slope of 
in this region reflects the ratio of the ~ 1 to 15 pm intensities of the EBL. 
The rise in beyond 10 TeV reflects the rise in the EBL towards the peak 
of the dust emission at ~ 100 — 200 pm. 

The energy dependence of r^.y will give rise to several breaks in the spec- 
trum of 7— ray sources that reflect the changes in the slope of the opacity. 
The first spectral break, ATcev occurs between 10 and 500 GeV. The second, 
Ar^ev around 1 TeV, and the third around 10 TeV. 

The first break has been used in most EBL studies to date with various 
assumptions on the intrinsic source spectra, and the second break has been 
most recently explored in the analysis of Orr et al. (2011). A review of studies 
that utilize the first and second break for constraining the EBL is presented 


9 




Figure 4: Left panel: Calculated EBL intensity versus wavelength at z = 0 ; Right 
panel: The 7 — ray opacity versus energy for sources at different redshift (see labels). The 
figure illustrates the correlation between the changes in the slope of the EBL intensity 
with those in Model calculations by Finke et al. (2010). Details in §2.5 of the text. 

in §6. A break at ~ 10 TeV has yet to be discovered. Starburst galaxies, 
which have a hard 7— ray spectrum, are the most promising subject for such 
analysis. 

3. The Types and Spectra of Extragalactic GeV /TeV Sources 

Determination of the EBL intensity from GeV-TeV 7— ray observations 
requires knowledge of the intrinsic spectrum of the sources. Here we list 
the different sources, their spectral characteristics, the different proposed 
mechanisms for their 7— ray production, and the physical limits on their 
spectral energy distribution at very high energies. 

3.1. The Types of Extragalactic GeV/TeV Sources 

The currently available GeV-TeV 7— ray sources that are being used to 
derive limits and constraints on the EBL are listed in Table 2. They in- 
clude the accretion-powered relativistic jets of active galactic nuclei (AGNs), 
namely BL Lacertae objects (BL Lacs), flat spectrum radio quasars (FSRQs) 
and a few radio galaxies. The list is complemented by the recent detections 
of two nearby starburst (SB) galaxies. In contrast to AGNs, their 7— ray 
spectrum is generated by the cumulative effects of cosmic-ray acceleration 
in shocks generated by a large number of supernova remnants (Volk et ah, 
1996). 

Blazars : To date, the most numerous sources used in EBL studies are 
blazars. Historically, they have been divided into two sub-classes based on 
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their optical properties: FSRQs, characterized by strong emission lines; and 
BL Lacs, characterized by weak or lack of emission lines. Because of the 
weakness of their emission lines, the redshift determination of BL Lac blazars 
has proven difficult or even impossible in many cases. The status of blazars 
as bright GeV and TeV sources arises from the fact that their relativistic 
jets are closely aligned with the observers line of sight. Consequently, the 
luminosity of a 7— ray emission region moving relativistically along the jet 
axis in the direction of the observer is strongly beamed, enabling its detec- 
tion at cosmological distances. Occasional strong flaring activity renders the 
following BL Lacs: PKS2155 (Aharonian et ah, 2007b); Mrk 501 (Catanese 
et ah, 1997; Abdo et ah, 2011a), and Mrk 421, (Gaidos et ah, 1996; Acciari 
et ah, 2011a; Aleksic et ah, 2011) the brightest TeV sources; and the fol- 
lowing FSRQs: 3C 454.3 (Donnarumma et ah, 2009), and 3C 279, (Wehrle 
et ah, 1998) the brightest GeV emitters in the sky. The flaring has provided 
high quality 7— ray spectra and has led to their detection at redshifts as far 
as 2: ~ 0.5 at TeV energies with lACTs, and as far as 2: ~ 3.2 at ~ 10 GeV 
energies with Fermi. 

The combined GeV-TeV observations of blazars makes it possible to study 
their spectra over a larger range of redshifts, thereby enabling the studies of 
the EBL over a wider range of wavelengths. GeV photons interact mainly 
with UV/optical photons, whereas TeV photons probe mainly the near- to 
mid-IR region of the EBL. Since the intensity of the EBL is much lower at 
UV energies, the universe is transparent to 7— rays below 10 GeV, becoming 
essentially opaque for TeV sources at redshifts of z > 0.5. The Fermi Gamma- 
Ray Space Telescope provides important probes of the UV region of the EBL, 
and the GeV transparency can be used to test evolutionary models of the 
EBL to relatively large redshifts (;^ > 1). 

Radio galaxies: The jets in radio galaxies are signihcantly misaligned with 
respect to the observer’s viewing direction, and thereby provide no relativistic 
Doppler boosting. This limits the detection of radio galaxies with current 
generation 7— ray telescopes to the local group and the Perseus galaxy cluster. 
Deep 7— ray observations of radio galaxies with GTA combined with spatially 
resolved studies in the radio, optical and X-ray will play an important role in 
understanding the physics of relativistic jets. These observations are likely 
to yield spectra up to ~ 10s TeV which will provide useful constraints on 
the EBL in the mid- and far-IR wavelength regions (Acciari et ah, 2009c, for 
the VERITAS Gollaboration, the VLBA 43 GHz M 87 Monitoring Team, the 
H.E.S.S. Gollaboration, and the MAGIG Gollaboration). 
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With sufficiently high spatial resolution, the 7— rays produced by IC scat- 
tering of CMB and EBL photons off the relativistic electrons of the lobes of 
radio galaxies can be used to set limits on their energy density in the imme- 
diate vicinity of these objects (Georganopoulos et ah, 2008, see §6 below). 
Starburst galaxies; The detection of starburst galaxies M82 (Acciari et al., 
2009b; Abdo et al., 2009) and NGC 253 (Acero et al., 2009), potentially 
opened a new wavelength regime for studying the EBL. The 7— rays in star- 
burst galaxies are generated by cosmic rays that are accelerated by a large 
number of supernova remnants, giving rise to hard 7— ray spectra that ex- 
tend to energies of 10s of TeV. The ~ 10 TeV opacity to nearby starburst 
galaxies is quite small, and about unity at energies of ~ 50 - 100 TeV. Nearby 
starbursts are therefore important probes of the EBL at far-IR (~ 100 /rm) 
wavelengths that cannot be probed by other 7— ray sources because of the 
relative softness of their spectra compared to those of SB galaxies. 

3.2. The Spectra of Extragalactic GeV/TeV Sources 

Over a sufficiently small energy range the blazar spectrum can be charac- 
terized by a power law, dN/dE oc E”'", with different indices, TceV and Trev, 
at GeV and TeV energies, respectively. An important characteristic of the 
observed spectra is the presence of a break, dehned as ATcev = hcev ~ hTev, 
occurring between GeV and TeV energies, the exact location depending on 
the source’s redshift. A source with an intrinsic spectrum characterized by 
a single power law out to TeV energies will have a value of ATcev = 0. 
Without any intergalactic absorption this value will remain constant with 
redshift. 

The spectral index Tq^v is obtained from a power law £t to the ~ 1 — 
10 GeV region of the spectrum which is unaffected by EBL absorption. If 
the intrinsic blazar spectrum is an extension of this power law to energies of 
~ 1 TeV, then any spectral break {TtcV > ^Gev) in the observed spectrum 
can be regarded as evidence for EBL absorption. 

A spectral break analysis of the amount of EBL absorption provides there- 
fore a powerful method for studying the EBL. It is a differential method that 
replaces knowledge of the intrinsic blazar spectrum with weaker requirement, 
namely that the power law representing the intrinsic blazar spectrum at GeV 
energies can be extended to TeV energies as well. 

Table 2 lists the values of TceV and TteV and the redshifts for all GeV and 
TeV detected blazars. Almost all sources exhibit a spectral break (ATcev < 
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Figure 5: The difference between FceV, the spectral index at GeV (Fermi) energies, 
and FxeV, the energy spectral index in the TeV regime (H.E.S.S, MAGIC, VERITAS) 
is shown as a function of their redshift. Red squares (radio galaxies), red stars (starburst 
galaxies), empty circles (HBLs, high-frequency peaked BL Lacs), blue downward triangles 
(intermediate- frequency peaked BL Lacs), filled circles (LBLs, low-frequency peaked BL 
Lacs), red upward triangles (FSRQs, flat spectrum radio quasars) indicate the different 
types of 7 — ray sources. 

0) at energies between 10 GeV and the 1 TeV. Fignre 5 depicts the depen- 
dence of AToev on redshift. The fignre shows a clear trend of increasing 
lATcevI with redshift, strongly snggesting that the break is the conseqnence 
of the attennation of the sonrce spectrnm by the EBL. As the optical depth 
increases with redshift, the observed 7 — ray spectrnm becomes softer, the 
position of the break moves to lower energies, and AToev becomes more neg- 
ative. The detailed redshift dependence of ATcev reflects the evolntion of 
the spectrum and proper photon number density of the EBL with redshift 
(see §7). 

Figure 5 also shows that there is significant scatter in AFcev at any given 
redshift. This suggests that some sources have considerable intrinsic harden- 
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ing/ softening in their spectra. Indeed, observations show that sources exhibit 
a wide range of spectral trends, i.e., spectral softening (3C 279, PKS 1510-08) 
and spectral hardening (lES 0502-1-675) in the GeV regime (Abdo, A. A. et 
ah, 2010), while the spectra of other blazars (Mrk 421, Abdo et ah, 2011b), 
(Mrk 501, Abdo et ah, 2011a) extend without any cutoff to 300 GeV (Abdo 
et ah, 2010b), or multi- TeV energies (Krennrich et ah, 2002), depending on 
the flaring state of the source. A lower limit in the scatter of AFoev at a 
given redshift may suggest the combined effects of intrinsic spectral softening 
of the source spectrum and the effects of EBL absorption. An upper limit in 
the scatter may simply indicate that the break in the spectral index is only 
created by EBL attenuation. 

Disentangling intrinsic spectral softening from the effects caused by the 
EBL is complicated, and requires a clear understanding of the physical pro- 
cesses that cause the intrinsic softening or hardening (which is very rare) 
between the GeV and TeV energy regions of the source spectrum. 

Multiwavelength observations are required to solve this problem. Such 
observations will provide conclusive tests of non-thermal 7— ray emission 
models, and constitute an important step towards achieving the ultimate 
goal of using the radio-optical- X-ray and low energy 7— ray spectra as a pre- 
dictor for the TeV spectrum. The application of this approach to different 
source classes, i.e., blazars, radio galaxies and starburst galaxies could pro- 
vide additional redundancy to help constrain the EBL. 

The current constraints on the UV/optical to mid-IR regions of the EBL 
have come predominantly from studies of blazars, since the vast majority of 
the extragalactic GeV/TeV 7— ray sources are FSRQs and BL Lacs. AGN 
population studies with Fermi reveal only a small number of non-blazar 
AGNs (Ackermann et ah, 2011). Gonsequently, we will focus in the follow- 
ing sections on EBL limits that were derived from studies of blazar spectra, 
which have already provided important limits on the near- to mid-IR spectra 
of the EBL. 

3.3. Phenomenology of Blazar Spectra and Models for TeV ''y— ray Production 

The non-thermal emission spectra of blazars generally exhibit two emis- 
sion peaks in izF^, the power emitted per unit logarithmic photon energy 
[see Figure (6), (Abdo et ah, 2011b)]. The peak in the radio- UV-X-r ay wave- 
band is unequivocally attributed to synchrotron radiation that is produced 
by ultra-relativistic electrons. The second peak, located at X-rays or 7— ray 
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Figure 6: The SED of Mrk 421 depicts the synchrotron peak at X-ray energies, and the 
inverse Compton peak at TeV energies. The figure was taken from (Abdo et ah, 2011b, 
Figure 11). 


energies, is commonly believed to stem from soft photons that were upscat- 
tered by the inverse Compton process to X- and 7— ray energies by the very 
same electrons responsible for the synchrotron emission. This mechanism for 
creating the second peak is often referred to as the synchrotron-self- Compton 
(SSC) mechanism. If the second peak includes upscattered photons drawn 
from other ambient radiation helds the mechanism is called the external- 
Compton (EC) mechanism. 

Alternatively, the second emission peak can be produced by a hadronic jet 
containing a signihcant amount of energy in ultra-high-energy (UHE) protons 
that subsequently interact with soft photons via py — )■ 7 t°, tt^ interactions. 
The pions decay giving rise to 7— rays, electrons, muons, and neutrinos. The 
strong magnetic helds required to collimate the hadronic jets also lead to 
the generation of considerable synchrotron radiation by protons and charged 
leptons in the pair cascade. The electrons from the pair cascade contribute 
to the lower energy synchrotron peak, whereas the muons together with the 
secondary photons from neutral pion decay contribute to the higher energy 
7— ray peak. For a recent review of blazar models see Boettcher (2010) and 
Dermer (2012). 

BL Lacs are classihed loosely by their synchrotron peak position. They 
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are referred to as low-frequency peaked BL Lacs (LBLs) if their synchrotron 
peak is at 10^^ < i^peak < 10^"^ Hz; as intermediate-frequency peaked BL 
Lacs (IBLs) if 10^^ < z/peak < lO^^Hz; and as high-frequency peaked BL Lacs 
(HBLs) if z/peak > lO^^Hz (Nieppola et ah, 2006). The corresponding 7— ray 
peak follows a similar pattern, with the peak energy progressively shifting 
towards higher energy as the synchrotron peak shifts to higher frequencies. 
LBLs peak at a few GeV (similar to FSQRs), IBLs in the tens of GeV, and 
HBLs generally peak beyond 100 GeV. 

As the 7— ray peak is found at increasingly higher energy, its luminosity 
also decreases relative to that of the synchrotron peak, a trend often referred 
to as the Fossati blazar sequence, (Fossati et ah, 1998). As a result, the multi- 
wavelength spectra of FSRQs and LBLs exhibit a prominent and large 7— ray 
luminosity, whereas in IBLs and in particularly HBLs, the 7— ray peak is typ- 
ically dwarfed by the synchrotron emission. This makes both FSRQs and BL 
Lacs very useful and complementary for EBL studies. FSRQs are extremely 
bright in the GeV regime and despite their low peak energy, some are still de- 
tectable in the sub-TeV regime with GeV and TeV telescopes. They are thus 
becoming increasingly useful for EBL constraints in the UV/optical/near-IR. 
HBLs are at the other extreme, their 7— ray peak can extend well into the 
multi-TeV regime, while their 7— ray to synchrotron luminosity ratio is much 
smaller. However, during strong flares some HBLs have been found to show 
a 7— ray dominated spectral energy distribution. They are therefore useful 
for constraining the EBL at near- to mid-IR wavelengths. 

In general, HBLs can be well described by basic SSG models where the 
ultrarelativistic electrons and their target photons are closely linked via syn- 
chrotron radiation. IBLs are better described by external Gompton (EG) 
models, that include a strong ambient photon held external to the blazar jet, 
thereby providing additional target photons for IG scattering (Acciari et ah, 
2009e). FSRQs, such as 3G 279, are difficult to ht with either model and 
may require the addition of a hadronic component.” 

The detailed htting of blazar spectra with models requires extensive mul- 
tiwavelength monitoring including optical. X-ray, and 7— rays. Results from 
such multi- wavelength campaigns carry the potential to reduce the spread 
in the AFGev(z) relation in hgure 5. While numerous successful multiwave- 
length campaigns have been reported in the literature (Abdo et ah, 201 lb, a; 
Aleksic et ah, 2010a; Fossati et ah, 2008; Abramowski et ah, 2011), the 
unequivocal interpretation of the spectral energy distributions with a clear 
prediction for the intrinsic TeV spectrum has not been possible. 
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In spite of the fact that blazar spectra are complex, it is possible to set 
limits on the behavior of their GeV/TeV spectra that are based on funda- 
mental physical limits imposed by energy losses in the particle acceleration 
and radiation processes. The most prominent limitation is the maximum 
hardness of the 7— ray spectrum which provides an important constraint on 
the EBL. The energy spectrum of electrons produced in models of diffusive 
shock acceleration in blazar jets (Maikov & O’C Drury, 2001), strongly con- 
strains the hardness of the resulting 7— ray spectra produced in SSC and EC 
models, limiting their power law index to values larger than T ^ 1.5. 

Even spectra that obey this limit are difficult to produce at higher energies 
where the Klein-Nishina effect softens the energy spectra substantially. The 
detection of relatively hard TeV spectra of blazars with redshift 0.1 — 
0.2 therefore came as a surprise, since the absorption corrected spectra are 
already reaching the T 1.5 limit for the minimal EBL imposed by the IGL 
(Aharonian et ah, 2006a; Levenson & Wright, 2008; Krennrich et ah, 2008; 
Ackermann et ah, 2011; Abdo et ah, 2010a). 

To which extent a problem of hard TeV spectra persists depends on the- 
oretical scenarios explaining these hard spectra. A solution to the problem 
was proposed by Katarzyhski et ah (2006). In their model, a high low-energy 
cutoff in the electron distribution could give the appearance of a hard 7— ray 
spectrum for a given energy regime. Other ideas by Aharonian et ah (2008) 
show that 7 — 7 absorption in the source due to narrow band emission from 
the AGN could lead to unusually hard TeV spectra. The emission produced 
by proton synchrotron radiation (Aharonian et ah, 2007a; Zacharopoulou 
et ah, 2011) combined with internal absorption at lower energies has been 
shown to produce spectra that exceed the hardness achievable by DSA. Other 
explanations avoid substantial EBL absorption by introducing axion-like par- 
ticles that couple with photons in intergalactic magnetic helds thus reducing 
the 7— ray opacity of the universe substantially (Sanchez-Gonde et ah, 2009). 
Mechanisms that would substantially weaken EBL limits derived from 7— ray 
observations are discussed in Section 8. 

3.4. Unphysical Blazar Spectra 

In spite of the difficulties and the complexity associated with modeling the 
intrinsic 7— ray spectra of blazars, models generally predict a common feature 
of blazar spectra: they follow a power law with curvature and/or exponential 
cutoff and overall can be described by a concave shape (curvature in energy 
flux is downward rather than upward). This means that any absorption 
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corrected gamma-ray spectrum showing an exponential rise cannot represent 
a physical source spectrum, and must have its origin in an over-correction 
for EBL absorption^. 

The relation between the intrinsic, (dN/dE)int, and observed, (dN/dE)obs, 
blazar spectrum is given by: 


dN 

He 


int 


dN 

He 


y(E^,z) 


obs 


( 11 ) 


The equation illustrates that an overestimate of the 7 — 7 opacity will lead 
to an exponential rise in the inferred intrinsic spectrum of the blazar. Such 
exponential rise is unphysical. It runs contrary to our basic understanding 
of blazar models, and is absent in the 7 — ray spectra of blazars at energies 
below 10 GeV (Abdo et ah, 2010a). This behavior can therefore be used to 
exclude EBL models with optical depths that will result in an exponential 
rise in the corrected blazar spectrum (Guy et al., 2000; Dwek & Krennrich, 
2005). 

Similarly, EBL scenarios can lead to an absorption correction for which 
the reconstructed intrinsic spectra follow a power law with an extremely hard 
spectral slope. Additional constraints can be derived from the spectral slope 
itself, however these are model dependent. EBL models leading to slopes of 
T < 1.5 for the intrinsic spectrum, have been rejected by Aharonian et al. 
(2006a) on the basis of the diffuse shock acceleration model. Gaveats to this 
approach have been extensively discussed in the literature (Stecker et ah, 
2007; Katarzyhski et al., 2006; Bottcher et al., 2008; Aharonian et al., 2008; 
Krennrich et al., 2008; Lefa et al., 2011; Zacharopoulou et al., 2011), and 
while some extreme blazars may exhibit harder spectra with T < 1.5, most 
GeV energy spectra of blazars obey the T = 1.5 limit, with few exceptions 
(Ackermann et ah, 2011). 

These general constraints on the hardness of the intrinsic blazar spectra 
have been used to derive constraints on the EBL by using few individual 
objects (Dwek & Krennrich, 2005; Aharonian et al., 2006a), and by using 


^This assumes that no other external effects such as pileups due to pair halos are 
playing a major role in the modification of blazar spectra. The production of pair halos 
would require a magnetic field intensity between 10”^ and 10“^^ G, sufficiently large so 
that electrons are isotropized, and sufficiently small so that IC losses exceed synchrotron 
losses. Instead of a power law with curvature term, a power law with exponential cutoff 
or a broken power law is also consistent with the above statement. 
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samples of blazars (Mazin & Raue, 2007; Orr et al., 2011). 


4. The Extragalactic Background Light I: Measurements 
and Limits 


4-1- Spectral Measurements and Limits 

The specific intensity of the EBL is usually presented in units of nW m~^ sr 
The conversion between these units and others that are sometimes used in 
the literature is given by: 


z/ J,^(A) [nW m ^ sr ^ 


3000 
9.85 X 10-3 


\{pm 


h{\) [MJy sr-3] 

Iy{\) [mJy deg-^] 


( 12 ) 


The EBL intensity can be determined in several ways. The hrst con- 
sists of direct measurements, a method that poses considerable technical and 
astronomical challenges. Technically, it requires the absolute calibration of 
the instruments, and the understanding and removal of all measurements 
uncertainties. Astronomically, it requires the removal of strong foreground 
emission from interplanetary dust particles (the zodiacal light, ZL) and from 
stellar and interstellar emission components in the Milky Way. A thorough 
review of the challenges in determining the EBL was presented by Hauser & 
Dwek (2001). 

A strict lower limit to the EBL intensity can be obtained by adding 
the light emitted by resolved galaxies. In principle, the integrated galaxy 
light (IGL) can converge to the total intensity of the EBL. A necessary 
condition for convergence is that the spectral index a of the differential galaxy 
number count versus flux S', dN/dS ~ S~°‘, becomes smaller than 2 at lower 
fluxes, so that the total integrated intensity, J {dN/dS) dS, is hnite. At 
short wavelengths the intensity of the IGL is limited by the sensitivity of 
the survey. However, even in deep surveys the convergence of the IGL does 
not ensure the measurement of the total EBL intensity, since the low surface 
brightness regions of galaxies may be missed in standard aperture photometry 
(Bernstein et ah, 2002; Levenson & Wright, 2008). Furthermore, a truly 
diffuse background will always remain undetected in such surveys. 

At longer wavelengths and large beam sizes, unresolved galaxies become 
a source of confusion, limiting the depth of the survey. Below a certain 
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flux (the confusion limiting flux) individual sources become indistinguishable 
from the variation in the sky brightness caused by statistical fluctuations 
in the number of faint resolved or unresolved sources (Dole et ah, 2004, 
and references therein). These limitations can be partially circumvented 
by stacking analysis. Stacking of astronomical images of sources detected 
at one wavelength enhances their signal relative to the random background 
fluctuations at some other wavelength (e.g. Dole et ah, 2006). The integrated 
light obtained by this method is thus closer to the EBL intensity than that 
obtained by integration down to the confusion limit. 

Finally, given a model for dN/dS, one can extrapolate the differential 
source count to very faint fluxes, and evaluate the sensitivity of the integrated 
intensity to the lower flux limit and functional shape of the extrapolation. 

Tables 3-5 list measurements and limits on the EBL intensity derived by 
the different methods described above with the different satellites, balloons, 
and ground observatories. Select measurements were used to dehne the gray 
area in Figure 7. Absolute measurements and their la uncertainties were 
used to dehne the upper limits on the EBL. The integrated light from resolved 
galaxies and their la uncertainty was used to dehne the lower limit on the 
EBL. Lower limits derived from stacking analysis were used when available. 
At 140 fim the DIRBE detection with the FIRAS calibration (Hauser et ah, 
1998), and at longer wavelengths the FIRAS detections by (Fixsen et ah, 
1998) were used to dehne the limits on the EBL. The measurements used 
to dehne the upper and lower limits on the EBL are shown as bold entries 
in Tables 3-5. The hgure shows that the EBL is poorly determined in the 
~ 5 — 60 /rm wavelength region, where the foreground emission from the 
interplanetary dust cloud is strongest (Kelsall et ah, 1998). 

4-2. Integral Constraints on the EBL Intensity 

The total EBL intensity per unit solid angle is given by: 

where C{z) is the luminosity density in a comoving volume element at redshift 
2 :, and 1^1 is given by eq. (3). 

The comoving luminosity density is dominated by the radiative output 
from stars. On a galactic scale, an AGN can dominate the optical to IR 
output of a galaxy, however, on a global scale AGN make only a small con- 
tribution to the total energy releases in the universe. AGN make up most of 
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dz 


dz 
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A(yU,m) 

Figure 7: Limits on the EBL intensity. Lower limits (blue squares) are determined by the 
intensity of the IGL. Upper limits (red circles) are determined by absolute measurements 
of the EBL. The data used in the figure are listed in Tables 3-5 in bold. The shaded area 
depicts the range of the allowed EBL intensity as determined by UV to sub millimeter 
observations. 

the X-ray background (Mushotzky et ah, 2000; Draper & Ballantyne, 2009), 
and a significant fraction of the cosmic radio background (Dwek & Barker, 
2002). However, they make only a small contribution to the total IR back- 
ground. Accretion onto a central black hole releases about 10% of the rest 
mass energy of the accreted matter, significantly more than the 0.7% re- 
leased in nuclear processes. However, the mass locked up in BH is only 
~ 0.6% of that in stellar objects. Rest frame color-color diagrams generated 
from Spitzer IRAC and MIPS observations covering the 3.6, 4.5, 5.8, 8.0, 
24, and 70 pm bands, of radio-detected submillimeter-selected galaxies with 
spectroscopic redshifts show that AGN constitute a small fraction, between 
13 and 19%, of the sample dominating its mid-IR spectrum (Hainline et ah, 
2009). So with these limits in mind, the total comoving luminosity density is 
a direct measure of the cosmic star formation rate (CSFR) at a given redshift. 

Figure 8 presents the redshift dependence of the CSFR as determined 
from UV-optical emission lines, [O II], [O HI], Lya, Ha, H/3, and from mid- 
IR, submillimeter, and radio observations (Madau et ah, 1996; Hopkins & 
Beacom, 2006; Michalowski et ah, 2010). The data in the figure were taken 
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redshift z 


Figure 8: A compilation of the cosmic star formation rate as inferred from UV, Hq,, mid-IR, 
submillimeter, radio, and Lya observations. Lower limits were excluded from the figure. 
The red curves are analytical approximations representing the upper and lower limits to 
the CSFR. 

from the compilation of Michalowski et al. (2010), who used the standard 
ACDM parameters to calculate volume densities and the Kennicutt (1998) 
law with a Salpeter stellar initial mass function (IMF) to convert luminosity 
densities to star formation rates. The two red lines represent a broken power 
law approximation to the upper and lower limits lines of the observations. For 
consistency with the observational determination of the CSFR, the CSFR was 
converted to a bolometric intensity using the Salpeter IMF and a starburst 
age of 100 Myr: Lboi = 7.5 x 10® T Mq yr“^ (Dwek et ah, 2011). The 
integrated intensity of the CSFR is then bounded by: 

Iebl = 21 — 66 nW m“^ sr“^ (14) 

Table 6 presents the current limits on the EBL intensity, separated into 
its stellar and dust components, and compares them to this predicted by the 
various EBL models presented in §5. 

4-3. Constraints on the EBL from Fluctuation Measurements 

Most of the EBL is generated by discrete galactic or primordial stellar 
sources. Fluctuations in their number and their clustering properties will give 
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rise to spatial fluctuations in the EBL intensity. Studies of the optical region 
of the EBL via the fluctuations method were first conducted by (Shectman, 
1973, 1974), and of the near-IR region using the C'OBE/DIRBE data by 
Kashlinsky et ah (1996). The fluctuations do not provide a direct measure- 
ment of the EBL intensity. The derivation of the EBL intensity from these 
measurements will require detailed knowledge of the galaxy source counts, 
their luminosity function, and their clustering properties as a function of 
redshift. However, spatial fluctuations in the EBL provide a different means 
of setting limits on its intensity. Fluctuation measurements of the EBL do 
not require absolute measurements of its intensity, since the removal of fore- 
ground emission components is done on the basis of their distinct spatial 
properties rather than their absolute intensities. At IR wavelengths, these 
fluctuations can have a spatial signal that is distinctly different from that 
generated by the interplanetary dust cloud or by interstellar dust. After the 
removal of all known resolved sources, any residual fluctuations will measure 
the EBL contribution from sources that may represent a yet unknown or 
unresolved population of stars or galaxies. Recent fluctuation measurements 
in the Spitzer 3.6 and 4.5 pm bands (Kashlinsky et ah, 2012) show that they 
are in excess of those that can be attributed to known galaxy populations 
(Helgason et ah, 2012), suggesting an origin in a very faint, yet unknown, 
population of highly clustered sources. More details and references on the 
use of fluctuations analysis to constrain the EBL intensity can be found in 
Hauser & Dwek (2001) and (Kashlinsky, 2005). 

EBL fluctuations at far-IR wavelengths were detected at 170 pm (Lagache 
& Puget, 2000, ISO data), at 160 pm (Lagache et ah, 2007, Spitzer data) and 
by (Shang et ah, 2012), using the Planck data. Fluctuation measurements 
and EBL colors have been used by Penin et ah (2012) to derive limits on the 
EBL at 100 and 160 pm (see Tables 4 and 5). 

5. The Extragalactic background Light II: Models 

The EBL intensity only provides an integral constraint on all the radiative 
energy releases over cosmic time. In a dust-free universe, it represents the 
intrinsic stellar or AGN spectra. The EBL intensity and spectral shape de- 
pends then on the star formation history, the stellar initial mass function, the 
evolution of metallicity, the energy released by AGN, and the relative impor- 
tance of energy releases by nuclear and gravitational processes. In a dusty 
universe, the comoving luminosity density, C{z), and total EBL intensity. 
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Iebl remain unchanged, however the energy is redistributed by absorption 
and reemission processes over a large spectral range. The resulting spectrum 
depends on many factors, ranging from the size distribution, composition, 
and optical properties of the dust grains, the evolution of their abundance 
and properties over time, and on the morphology of the galaxy which deter- 
mines the spatial distribution of the dust relative to the radiative sources. 

Several distinct approaches have been used to model the intensity and 
spectral distribution of the EBL at ;^ = 0. They all represent different ap- 
proaches for calculating the evolution of Cy{\,z) with redshift [see eq. (2)]. 
Backward evolution models start from the local determination of C^{X,z = 
0), evolving it with redshift using observed galaxy number counts at different 
wavelengths. Forward evolution models use the CSFR to determine C^{\, z) 
as a function of redshift, and population synthesis and radiative transfer 
models to determine the distribution of the energy over wavelengths. Cosmic 
chemical evolution models are similar to previous models, except that their 
system is the universe as a whole. Semi-analytic models calculate £,y(A, z) 
by including the appropriate physical processes in a more general model for 
the formation and evolution of structure in the universe. In the following we 
discuss these different models in somewhat more detail. For a more exten- 
sive discussion and description we refer the reader to the review by Hauser 
& Dwek (2001). 

5.1. Backward Evolution (BE) Models 

BF models start with the construction of a library of the galactic spec- 
tral energy distributions (SFDs) representing those of galaxies in the lo- 
cal universe and evolve them back in time in order to £t observed number 
counts. The FBL intensity and spectral shape is used as an integral con- 
straints on their evolution. The galaxies in such library should represent the 
range of observed galactic morphologies (spiral elliptical, irregular) and ac- 
tivities (AGN, normal, starburst, mergers) in the local universe. The SFDs 
of the galaxies comprising such library should also satisfy observed statistical 
properties of the ensemble of local galaxies, such as: the trend of increas- 
ing S {60 fim)/ S {100 fim) /) and decreasing S {12 jj,m) / S {25 jam) /) flux ratios 
with increasing IR luminosities ((Soifer & Neugebauer, 1991); and the num- 
ber density of galaxies in the L + dL luminosity interval, represented by the 
luminosity function (LF), ^{L)dL. 

Many functional forms have been adopted to characterize the local LF. 
The most commonly one used at optical and near-IR wavelengths is the 
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Schechter (1976) LF. At wavelengths above ~ 10 /im the galaxies’ SED is 
dominated by thermal emission from dust, and their LF cannot be adequately 
represented by the Schechter LF. Several distinct functional forms have been 
used to characterize the LF in the different IR wavebands. They are usually 
characterized by three parameters: a normalization parameter, $*, a char- 
acteristic luminosity, L*, that determined the transition point between the 
low and high luminosity behavior of the L; and a power law index, a that 
determines the behavior of the LF at low luminosities. A list of references to 
the functional forms derived from the IRAS survey can be found in Hauser & 
Dwek (2001). Of those, the parameters of the Saunders et ah (1990) LF have 
been recently updated to ht the differential 24 /rm number counts obtained 
by deep Spitzer surveys (Rodighiero et al., 2010). 

If neither the galaxies’ SED nor their comoving number density evolved 
with time, the spectral luminosity density, would be independent 

of redshift. However, the recent deep surveys with the Spitzer and Herschel 
satellites show strong evolution in number counts, compared to predictions 
made with no evolution models (Lagache et al., 2005; Rodighiero et al., 2010). 
Evolution in the LF or, equivalently, the spectral luminosity density can be 
inferred directly from observations if the redshift of the sources is known, 
and their number counts are complete (e.g. Rodighiero et ah, 2010; Dunne 
et ah, 2000). Alternatively, evolution in the LF can be introduced by adding 
a redshift dependence, usually of the form (1 -|- z^, to the basic parameters, 
$*, L*, and a, that characterize the LF. The value of 7 is then derived by 
htting model prediction to the observed galaxy number counts in a given 
waveband. Evolution in the LF can also be modeled by evolving the rela- 
tive number of the different type galaxies: quiescent star forming galaxies, 
starbursts, AGNs, and ellipticals with redshift (Rowan- Robinson, 2001, 2009; 
Lagache et al., 2003; Dominguez et ah, 2011), or by simply evolving the rel- 
ative number of ultraluminous infrared galaxies (ULIRGs), characterized by 
IR luminosities in excess of 1012 Lo, relative to the rest of the galaxy 
population [e.g. Ghary & Elbaz (2001)]. Having determined the evolution 
of the spectral luminosity density with redshift, the EBL is obtained by a 
simple integration of C^{\,z) over redshift. 

BE models are relatively simple, and their predictions can be easily com- 
pared to observed galaxy number counts, their magnitude-color and magni- 
tude number density relations, and their redshift distribution. They are only 
loosely constrained by physical processes. The predicted dependence of the 
comoving bolometric luminosity density with redshift should be consistent 
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with that inferred from limits and observations of the cosmic star formation 
rate. 

5.2. Forward Evolution (EE) Models 

FE models use the redshift dependence of the cosmic star formation rate 
(CSFR), inferred from a variety of wavebands and line observation (Madau 
et ah, 1996; Hopkins & Beacom, 2006), as a starting point in their calcu- 
lations. The CSFR is determined from the UV-optical line and continuum, 
and IR and radio emission in the various wavebands using statistically deter- 
mined conversion factors derived from galaxies in the local universe (Kenni- 
cutt, 1998; Haarsma et al., 2000). Determination of the CSFR is complicated 
by extinction effects at UV and optical wavelengths, and by the implicit as- 
sumption that the IR luminosity is powered by stars and representative of 
the total bolometric luminosity of the galaxies. Even if the total bolometric 
luminosity of a given galaxy is determined, its conversion to a star formation 
rate requires knowledge of the stellar IMF, a poorly determined quantity at 
high redshifts, and the duration of the starburst activity at each redshift. 
Once the CSFR is determined, FE models use population synthesis models 
such as PEGASE Fioc & Rocca-Volmerange (1997), Starburst99 Leitherer 
et al. (1999), Bruzual & Chariot (2003) or Maraston & Stromback (2011) to 
calculate the stellar bolometric and spectral luminosity density as a function 
of redshift. 

The most difficult part of this approach is determining the fraction of 
starlight that is absorbed by dust, and the spectrum of the reradiated IR 
emission. The SED of a galaxy can be determined with radiative transfer 
models, such as GRASIL (Silva et ah, 1998), DIRTY (Gordon et al., 2001), 
or DUSTY (Nenkova et al., 2000). Alternatively, one can use a parametric 
approach in which the fraction of UV-optical light absorbed by the dust, and 
the reradiated infrared spectrum are statistically determined from observa- 
tions. EBL spectra derived from FE evolution models were presented by 
Dwek et al. (1998), Razzaque et al. (2009), and Finke et al. (2010). 

Population synthesis models, combined with simple radiative transfer cal- 
culations, are useful for determining the UV to radio SED of individual galax- 
ies. However, the cosmological application of such models assumes that star 
formation is a monolithic process, in which in all galaxies star formation com- 
menced at the same redshift and evolved quiescently until the current epoch. 
The models do not allow for galaxy interactions, stochastic star formation 
histories associated with merger events, or any morphological evolution of 
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galaxies. Any discrepancies between model predictions and galaxy number 
counts must be introduced in an ad hoc fashion by evolving the stellar IMF, 
or by introducing a new population of galaxies at the appropriate redshifts. 

5.3. Cosmic Chemical Evolution (CCE) Models 

CCE models treat the universe as a closed system in which all galax- 
ies within a large comoving volume element are represented by their basic 
ingredients: stars, interstellar gas, metallicity, and radiation. Chemical evo- 
lution equations, analogous to those used to follow the chemical evolution 
of the Galaxy (e.g. Audouze & Tinsley, 1976; Tinsley, 1981; Pagel, 2001), 
are used to follow the evolution of the average stellar, gaseous, and radiative 
contents in each comoving volume in a self consistent manner. CCE mod- 
els were pioneered by Pei & Fall (1995), and most recently updated by Pei 
et al. (1999). Inputs parameters for their model are the mean rest frame 
UV luminosity density as a function of redshift, and the mass of the ISM 
gas as determined from H I column densities derived from studies of quasar 
absorption lines through damped Lyo; systems. The decrease in the ISM gas 
with redshift and the UV luminosity density were used to derive a solution 
for the evolution of the CSFR with redshift which is consistent with that de- 
termined from the extinction-corrected Ha, and with SCUBA 850 and ISO 
15 pm surveys. Similar to FE models, population synthesis models were 
then used to calculate the stellar SED at each redshift, and an LMC extinc- 
tion law was adopted to calculate the fraction of starlight absorbed by the 
dust. A power-law distribution in dust temperature was used to calculate 
the spectrum of the reradiated IR emission. The model reproduced various 
observational constraints, including the comoving rest-frame 0.44, 1.0, and 
2.2 pm spectral luminosity densities in the ~ 0 — 2 redshift interval, the 12, 
25, 60, and 100 pm local luminosity densities; and the mean abundance of 
metals in damped Lya systems in the ~ 0.4 — 3.5 redshift interval. 

5 . 4 . Semi- analytical (SA) Models 

SA models follow the formation and evolution of galaxies in a cold dark 
matter Lambda dominated (ACDM) universe using the cosmological parame- 
ters derived from the 5-year Wilkinson Mierowave Anisotropy Probe ( WMAP5) 
observations (Hinshaw et ah, 2009) as the initial conditions. SA models then 
follow the growth and merger of dark matter halos, and the emergence of 
galaxies which form as baryonic matter falls into the potential wells of these 
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halos. The fate of the infalling gas is determined by many different pro- 
cesses: the formation of stars in a multiphase interstellar medium, AGN and 
supernovae feedback processes that quench their formation, the evolution of 
the stellar radiation held, the heating and cooling of the interstellar medium 
and its chemical enrichment, the exchange of material with the intergalac- 
tic medium through infall and galactic winds, and the growth of the central 
black hole. A description of recent developments and references to previous 
work can be found in Somerville et ah (2011). Model prediction are compared 
to a basic set of observational constraints such as the observed characteris- 
tics of galaxies: their morphology, colors, and spectral energy distribution, 
and morphology; and their integrated cosmological properties: their number 
counts and luminosity function in different wavebands and redshifts, their 
mass function, the cosmic star formation rate, and the EBL generated by 
them. As in all EBL models, determination of the galaxies’ SED is com- 
plicated by the detailed microscopic and large scale parameters needed to 
calculate the amount of starlight that is absorbed by dust, and the spectrum 
of the reradiated emission. Recent SA models have combined the models for 
galaxy formation with radiative transfer models to determine the galaxies’ 
SED (Eontanot et ah, 2009; Eontanot & Somerville, 2011; Somerville et ah, 
2011; Younger & Hopkins, 2011). 

SA models are inherently complex, incorporating a large number of phys- 
ical processes, some poorly known, to derive galaxy properties. However, 
they are the most physically motivated models, and quite successful in re- 
producing a large number of observational constraints. 

5.5. Comparison of Model Predictions with Observations 

A detailed comparison of all model types with EBL limits and observa- 
tions was presented by Hauser & Dwek (2001). Here we will represent mostly 
the models that have been developed since then: BE models by Stecker & 
Scully (2006), Franceschini et al. (2008), and Dominguez et al. (2011); The 
EE model of Finke et al. (2010); and the SA model of Gilmore et al. (2011). 
Figure 9 compares the various models to the current limits and observations 
of the EBL. In general, all models, except for the BE models of Stecker et 
al. provide adequate hts to the EBL. 
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Figure 9: Models of the EBL are compared to observational limits on the EBL. 

6. EBL Constraints from 7— ray Observations of Blazars 

The attenuation of 7 — rays by the EBL can in principle be used to deter- 
mine the EBL intensity at wavelengths corresponding to the 7 — ray observa- 
tions. Neglecting the possible scattering or production of second generation 
7 — ray photons along the line of sight to the blazar, the intrinsic 7 — ray flux 
from the blazar, Fjnt(E^), can be related to the observed one, Fobs(E^) by: 

Eobs{E^) = Eint{E^) exp[-r.^.^(E^)] (15) 

where the optical depth, is given by eq. (9). Determination of the EBL 
assumes that all the attenuation is caused by interaction with the EBL, in- 
stead of photons in or around the vicinity of the blazar. Furthermore, it 
requires knowledge of the intrinsic blazar spectrum. Assuming that all the 
attenuation is attributed to the EBL, several upper limits have been derived 
on the EBL intensity by making various assumptions on the intrinsic blazar 
spectrum. The 7 — ray derived EBL limits are compared to those derived 
from UV to sub millimeter observations in Figure 10, and described below. 
Fixed power law: Early observations of blazars (Mrk 421 Punch et ah, 
1992) suggested that the ~ GeV-TeV spectrum could be approximated by 
a single power law. If so, then any deviations of the observations from the 
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extrapolated power law to higher energies should be attributed to EBL at- 
tenuation. Stecker & de Jager (1993) derived an upper limits at 1-5 pm 
of 10 nW sr“^, assuming that a straight power law extrapolation of 
the spectrum of Mrk 421 from GeV energies obtained from the EGRET 
(Lin et al., 1992) with an index of T = 1.96 ± 0.14 holds up to the TeV 
regime, where the index was measured by the Whipple collaboration to be 
r = 2.25 ± 0.19. Biller et al. (1995) included the statistical uncertainties 
of the GeV spectrum, and demonstrated vastly different extrapolations with 
signihcantly higher upper limits, thereby yielding conservative upper limit to 
the EBL in the mid-IR at 10pm. 



Figure 10: Limits on the EBL as determined from 7 — ray observations of blazars. Details 
in text. 

Strong flaring activity of Mrk 501 provided a well measured TeV spectrum 
from 0.2 - 24 TeV (Aharonian et ah, 1999; Samuelson et ah, 1998; Djannati- 
Atai et al., 1999). Stanev & Franceschini (1998) calculated limits by htting 
the energy spectrum of Mrk 501, with a range of possible absorption scenarios 
assuming intrinsic power law spectra and varying levels of EBL intensity by 
scaling the lower limits from galaxy counts. This provided strong EBL limits 
in the near- and mid-IR. Funk et al. (1998) followed a similar approach, 
except that they used an EBL model by MacMinn & Primack (1996) as the 
basis for scaling the EBL intensity, yielding similar results in the mid-IR. 
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Mannheim (1998) argued that the observed spectra would deviate from a 
power law, if the primary 7 — ray spectrum were substantially attenuated by 
the EBL. An upper limit in the mid-IR that is based on this hypothesis and 
an energy spectrum of Mrk 501 from the HEGRA collaboration is also shown 
in Figure 10. Later on, Vassiliev (2000) demonstrated that the absence of 
deviations from a power law does not preclude the presence of substantial 
absorption in the observed spectra. 

The obvious drawback of the method is that the assumption of a single 
power law for the intrinsic blazar spectrum does not hold true over a wide 
range in energy. Blazar spectra generally exhibit a concave 7 — ray peak over 
a sufficiently large energy range. Most blazar spectra measured by Fermi or 
by Cherenkov telescopes can be represented by a power law over the energy 
range covered by the instrument. 

Synchrotron Self-Compton (SSC) Spectrum; The SSC model is a pop- 
ular model explaining the existence of the two peaks in the blazar spectrum: 
the synchrotron peak at radio- UV-X-r ay energies and the inverse Compton 
(IC) peak at 7 — ray energies. The spectrum of the IC peak can be modeled 
using parameters that produce the synchrotron peak and the unabsorbed 
part {E < 10 GeV) of the IC spectrum (see review by Dermer, 2012). 

Such models for the intrinsic blazar spectrum have been used by Guy 
et ah (2000) to determine the intensity of the EBL in the 1-5 pm and 20-80 
pm wavelength region. They applied a multiwavelength £t to the X-ray and 
TeV data of Mrk 501 in the framework of a standard homogeneous SSC model 
to derive the level of absorption present in the TeV spectrum. As a result, 
they obtained an absolute upper limit on the EBL of 60 nW m“^ sr“^ and a 
most likely value of 20 nW m“^ sr“^ at 1 pm. They also pointed out that the 
lack of an absorption signature in the spectrum of Mrk 501, as suggested by 
the HEGRA telescopes, does not necessarily imply a lack of EBL absorption. 
They emphasized that in the transition region from the near-IR to the mid-IR 
EBL, the opacity could be nearly constant. This is a consequence of the large 
width of 7 — 7 cross section (see Figure 2 ). So when is convolved with the 
number density of background photons, any strong wavelength variations in 
the EBL are smoothed out. As a result, the observed TeV spectrum at 1 - 10 
TeV would corresponds to the intrinsic blazar spectrum since the observed 
spectrum is now described by (dN / dE)int x e^ with r a slowly varying function 
of energy. 

The drawback of using 7 — ray emission models to constrain the EBL is the 
uncertainty in the many parameters that determine the IC spectrum. Fur- 
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thermore, while HBLs generally can be well fit by SSC models, IBLs require 
the inclusion of additional ambient radiation helds that make a contribution 
to the 7— ray IC component. 

Additional complications arise from the fact that basic one-zone SSC 
models are not applicable for sources exhibiting ’’orphan flares”, where only 
the TeV flux is enhanced while the synchrotron emission remains unchanged 
(Krawczynski et ah, 2004). Finally, the biggest challenge for the SSC/multi- 
wavelength approach to constraining the EBL is to get simultaneous mea- 
surements for large sets of blazars. 

The r > 1.5 limit on the hardness of the blazar spectrum; A more re- 
laxed assumption on the intrinsic blazar spectrum is that it cannot pro- 
duce too many hard photons, so that the 7— ray spectrum, expressed as 
E'^dN/dE ~ E~^ cannot be flatter than one with F = 1.5 (Maikov & O’C 
Drury, 2001). In the spirit of this limit to the spectral index, Renault et al. 
(2001) explored a range of EBL scenarios based on measurements with the 
minimal assumptions that the intrinsic power of Mrk 501 is concave, effec- 
tively requiring a decreasing energy flux distribution above 4 TeV (F > 2.0). 
They derived an upper limit of 5 nW/m^/sr at 10 /im. 

The strict assumption of F > 1.5 , was used by Aharonian et al. (2006a) to 
derive upper limits on the 1-5 /rm on the EBL which are close to the lower 
limits determined by the IGL, suggesting that the EBL has been largely 
resolved at these wavelengths. A comprehensive study by Mazin & Raue 
(2007) is based on eleven blazars over a redshift range from 0.03 - 0.18, and 
explores a large number (8 million) of hypothetical EBL scenarios to set 
upper limits on the EBL, again with the requirement that the source spectra 
cannot be harder than F = 1.5orF = 2/3. The lower value arises from the 
extreme scenario of a mono-energetic energy distribution of ultra-relativistic 
electrons in which the resulting IC 7— ray spectrum could be as hard as 
7 = 2/3, leading to two conditional upper limits. The first condition yielded 
limits that are slightly above that of Aharonian et al. (2006a). The second, 
more relaxed condition, yielded limits that were higher by about 30%. 

The theoretical validity of a strict hardness limit of F > 1.5 has been 
discussed by a number of authors, with no unanimous verdict (Katarzyhski 
et ah, 2006; Stecker et ah, 2007; Bottcher et ah, 2008; Aharonian et ah, 
2008; Lefa et ah, 2011; Zacharopoulou et ah, 2011). Observational evidence, 
e.g. Levenson & Wright (2008), have provided lower EBL limits from galaxy 
counts that are higher than previous ones derived by Madau & Pozzetti 
(2000). If these new limits are correct, they imply 7— ray spectra that are 
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slightly harder than F = 1.5 (Krennrich et ah, 2008). 

Unphysical exponential rise of the blazar spectrum: Less model de- 
pendent, and therefore more robust limits on the hardness of the intrinsic 
blazar spectra arise from the notion that an exponential increase of their 
luminosity with energy is unphysical. All current blazar models produce 
a concave spectrum, rendering intrinsic blazars spectra with an exponential 
rise in energy flux theoretically unfeasible. The paradigm of concave intrinsic 
energy spectra was used by Dwek & Krennrich (2005) to reject many differ- 
ent realization of the EBL. Furthermore, Dwek et al. (2005b) ruled out the 
extragalactic origin of the near-IR sky brightness observed by Matsumoto 
et al. (2005), since it would lead to an exponential rise in the spectrum of the 
blazar PKS 2155-304, which is ruled out by observations (Aharonian et al., 
2005b). An EBL spectrum close to the IGL limits yielded a blazar spectrum 
consistent with the SSC model, suggesting that the EBL was mostly resolved 
at near-IR wavelengths (Dwek et al., 2005b, a). 

Spectral break analysis due to EBL spectrum: Orr et al. (2011) de- 
veloped a novel approach to set limits on the EBL intensity by studying the 
effects of the spectral shape of the EBL on the sub- to multi-TeV spectra 
of blazars. The 7 — ray opacity around 1 TeV is sensitive to the EBL inten- 
sity at ~ 1 //m (see Fig. 3). Its subsequent energy dependence hinges on 
the rate at which the stellar UV, optical, and near-IR emission decreases 
towards mid-IR wavelengths. The near- to mid-IR intensity ratio of the EBL 
determines the relative ~ 1 to 10 TeV opacity. This is illustrated in Figure 5 
which shows the energy dependence of A large near- to mid-IR ratio in 
the EBL intensity would cause the 7 — ray opacity to be relatively flat in the 
~ 1 — 10 TeV region, resulting in a hard 7 — ray spectrum. Conversely, a low 
near- to mid-IR ratio would result in an increase in the ~ 10 TeV opacity 
relative to that at ~ 1 TeV, resulting in a softer blazar spectrum. 

The spectral shape of the EBL in the ~ 1 — 15 /rm region can be related 
to the break AFrev = F(< 1 TeV) — F(> 1 TeV) in the blazar spectra. 
A study of 12 blazars spanning a redshift range from 0.03 to 0.186, showed 
a trend of increasing lAF-pevI with redshift with a statistical significance of 
3.6(7 (Orr et al., 2011). This strongly suggests that the trend is caused by 
EBL absorption, providing strong constraint on the near-IR to mid-IR EBL 
intensity ratio. Combined with 7 — ray-derived upper limits on the EBL (Orr 
et al., 2011) derived correlated constraints between the near-IR and mid-IR 
intensity ratio of the EBL and the EBL intensity at mid-IR wavelengths. 

7 — ray inverse Compton emission from radio lobes: A new method for 
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determining the local EBL intensity was described by (Georganopoulos et ah, 
2008). The method relies on the detection of 7 — rays produced by IC scat- 
tering of CMB and EBL photons off the relativistic electrons of the lobes of 
radio galaxies. Since the lobes have to be clearly resolved with the 7 — ray 
telescope, the method is presently limited to the nearby radio galaxy Eor- 
nax A. With the normalization and maximum electron energy in the lobes 
determined by their synchrotron spectra (measured by WMAP), the 7 — ray 
spectrum of the lobes is determined by the intensity of the radiation field 
in the radio lobes, which is dominated by the CMB and EBL. The result- 
ing 7 — ray spectrum comprises distinct contributions by the CMB and EBL 
photons: The CMB contribution peaks at energies ~ 40 MeV, while that of 
the EBL appears as an excess above steeply dropping CMB contribution at 
higher energies. This excess emission is relatively flat, and extends to ener- 
gies of ~ 50 CeV. Its magnitude provides a direct measurement of the EBL 
intensity at mid- and far-IR wavelengths. 

Combined 7 — ray limits: Eigure 11 shows select limits on the EBL de- 
rived from TeV observations. All limits were adjusted to a common Hubble 
constant of 70 km s“^ Mpc“^. Strict convergence of the different limits is 
not expected, since each limit was derived from different assumptions on the 
intrinsic blazar spectra and the EBL spectrum. The bold horizontal lines and 
the Orr et al. (2011) limit represent the most recent limits on the EBL. The 
1 /im and lO/rm upper limits derived in these studies agree within 20%. It is 
important to emphasize that different methods were used, i.e., the T > 1.5 
limit on the hardness of the blazar spectrum and the spectral break analysis 
due to EBL spectrum, yet they reached similar conclusions. 

In general, much improved upper limits are now available from the ob- 
servations with the new generation of atmospheric Cherenkov telescopes 
(H.E.S.S., MACIC, VERITAS). Both the near-IR and the mid-IR intensity 
levels are now constrained to much lower values than was possible with early 
results from the first few TeV blazars, and make only minimal assumptions on 
the intrinsic blazar spectra. Extreme EBL scenarios, with 70 nW m~^ sr“^at 
1.5 /im, as suggested by Matsumoto et al. (2005) are clearly ruled out. Min- 
imal assumptions on the blazar spectra such as the absence of exponential 
rises or applying the T > 1.5 limit, reject such high EBL intensities in the 
near-IR (Dwek et ah, 2005b; Aharonian et ah, 2006a). 

Mazin & Raue (2007) provide upper limits that are comparable to the 
former two in the near-IR, while in the mid-IR their upper limits are much 
higher. Eor this work the spectrum of lES 0229-1-200, a blazar with a redshift 
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Figure 11: The EBL constraints as derived from 7 — ray observations to the EBL are shown 
for 1 fjLm (left) and for 10 /rm (right). The horizontal line represents the intensity of the 
integrated galaxy light (IGL), which provides a strict lower limit on the EBL intensity. 
The most recent limits are represented by bold lines. The shaded area represents the Icr 
uncertainty in the IGL intensity. References to the 7 — ray limits are: (0) Stecker & de 
Jager (1993); (1) Stanev & Eranceschini (1998); (2) Biller et al. (1998); (3) Guy et al. 
(2000); (4) Dwek et al. (2005b); (5) Schroedter et al. (2005); ( 6 ) Aharonian et al. (2006a); 
(7) Mazin & Raue (2007); ( 8 ) Finke & Razzaque (2009); (9) Orr et al. (2011); (10) Funk 
et al. (1998); (11) Mannheim (1998); (12) Dwek & Slavin (1994); (13) Renault et al. 
(2001); (14) Aharonian et al. (2007a). 


of z=0.129 and a spectrum up to 10 TeV was not available, while it is the basis 
for strong mid-IR limits in Aharonian et al. (2007b) and Orr et al. (2011). 
This emphasizes the importance of extending the TeV energy spectra of 
distant (z> 0.1) blazars into the multi-TeV regime, where 7 — rays reach their 
maximum cross-section with photons in the mid-IR. Furthermore, energy 
spectra between 100 GeV and 10 TeV are sensitive to the spectral shape of 
the EBL, thereby linking the upper limits in the near- and mid-IR; a given 
mid-IR intensity level combined with the EBL spectral shape limits the range 
of near-IR EBL intensities, i.e., by excluding EBL intensities that are either 
too high or too low in the near-IR based on the spectral shape constraint. 

There are now in the order of 30 extragalactic objects with redshifts up 
to ^ ~ 0.5 available for constraining the EBL. Limits derived from the recent 
generation of experiments (H.E.S.S., MAGIC, VERITAS) suggest a low mid- 
IR, such as the results from Aharonian et al. (2007b) and Orr et al. (2011). 

In summary, it should also be noted that all of the most recent TeV 
constraints are well within the boundaries set by direct measurements and 
their uncertainties (see shaded area in Figure 10). However, given the debate 
about some of the assumptions about the blazar spectra, namely, the F > 1.5 
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limit on the hardness of the blazar spectrum, claims that the EBL has been 
resolved are premature. 


7. The 7 — ray Opacity of the Universe 


With blazars being detected at increasingly large redshifts, it becomes 
possible to use them to discriminate between different EBL models. The 
7— ray opacity to a blazar at redshift is given by: 




-log 


F(E^)int 


(16) 


The opacity can be determined from models of the EBL if its evolution with 
redshift is known and, independently, from 7— ray observations if the intrin- 
sic blazar spectrum is known. Concordance between these two independent 
determinations of Tjj can serve as a test for the validity of the underlying 
assumptions in each method. 

Eigures 12 and 13 depict the evolution of the comoving intensity of the 
EBL, the corresponding evolution of the proper number density of back- 
ground photons, the optical depth to blazars at various redshift, and the 
corresponding attenuation factor. Results are plotted for the BE evolution 
model of Eranceschini et al. (2008) and the BE evolution model of Dominguez 
et al. (2011). 

Determining the 7— ray opacity from observations requires knowledge of 
the intrinsic blazar spectrum. Differences between the observed and expected 
flux at a given energy would then be simply attributed to EBL attenu- 
ation. Eigures 12 and 13 show that the sharp drop of the EBL intensity at 
UV and shorter wavelengths renders the universe almost transparent to GeV 
photons. Consequently, the observed ~ 1 — 50 GeV spectrum is very likely 
the intrinsic blazar spectrum. So instead of assuming a theoretical limit on 
the spectral index, one can use the GeV - 10s of GeV energy spectral slope 
from Fermi data as a proxy for the intrinsic spectra at TeV energies. 

Assuming that this power law can be extrapolated from GeV to TeV 
energies, one can derive the TeV optical depth to the observed blazar. This 
approach was used by Georganopoulos et al. (2010) and in method 1 in 
Orr et al. (2011) to set hrm upper limits on EBL models using the GeV to 
TeV spectra of PKS 2155-304 {z = 0.116) and lES 1218-1-304 {z = 0.182). 
Assuming that the GeV spectrum is unattenuated by the EBL, Mankuzhiyil 
et al. (2010) used optical. X-ray and GeV data to model the TeV flux of 
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Figure 12: Basic EBL model results by Franceschini et al. (2008): Top left: The comoving 
EBL and CMB intensities versus wavelength for different redshifts; Top right: The proper 
number density of EBL and CMB photons versus energy for the same grid of redshifts 
as the previous panel; Bottom left: The 7 — ray opacity versus energy, for different 
redshifts; Bottom right: The amount of attenuation versus energy for the same grid of 
redshifts as the previous panel. The figure illustrates the change in the slope of at 
energies corresponding to the wavelength at which the slope of the EBL spectrum changes. 
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Figure 13: Same as Figure 12 for the Gilmore et al. (2011) model. 
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PKS2155-304 using a one- zone SSC model. Comparison of the model results 
with observations, they derived the TeV opacity to this blazar, and found it 
to be consistent with most EBL models. 

Figure 14 compares the dependence of the optical depth derived from 
EBL models (hatched curves) to that derived for select blazars: Mrk 501, 
lES 1218-1-304, and 3C 66A. Each hatched band spans the range of optical 
depths predicted by the EBL models of Franceschini et al. (2008), Finke et al. 
(2010), Dominguez et al. (2011), and Gilmore et al. (2011). The colored 
dots represent the optical depths derived from the 7— ray observations of the 
three blazars. The intrinsic blazar spectrum was assumed to be a power 
law determined by the observed flux at 1 GeV and the spectral index, TceV- 
The observed flux in the TeV range was assumed to be a power law with 
a spectral index TtcV (see Table 2). The 7— ray opacity in the TeV range 
was then derived from eq. (9). The band of opacities for each blazar was 
obtained by performing 100 Monte Carlo simulations of the intrinsic and 
observed spectra using the uncertainties in the spectral indices and 7— ray 
energies into account. 

The hgure shows that the 7— ray derived optical depths of Mrk 501 and 
lES 1218-304 are in general agreement with model prediction. The discrep- 
ancy between the the EBL and the 7— ray derived optical depth for 3C 66A 
is typical of most blazars listed in Table 2. The convergence between ob- 
servational limits on the EBL and models suggests that the origin of the 
discrepancy can be mostly attributed to our still incomplete knowledge of 
the intrinsic spectra of blazars. 

8. Is the Blazar Spectrum Determined in the Source? 

In deriving upper limits on the EBL from TeV observations, it was tacitly 
assumed that blazar spectra are produced in the sources and attenuated 
en route to Earth. The detection of medium redshift blazars at ~ 0.2 
with hard 7— ray spectra has revived alternate models for the origin of their 
spectra. Since hard photons from these redshifts are easily attenuated, one 
of these models proposes that the observed blazar spectra contains secondary 
photons that are produced close to the observer, and therefore more likely 
to survive. 

High resolution optical spectra of the blazars PKS 0447-439 and PMN J0630- 
24 obtained with the CTIO and NTT observatories show a clear absorp- 
tion line around 6280 A which, when attributed to the Mg II 2795.5 A and 
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Figure 14: Limits on the optical depth at z = 0.1, 0.3, 0.5, and 1.0 as determined by the 
EBL models are compared to observationally determined optical depths for select blazars. 
The optical depth to 3C 66A is still uncertain. See §7 for more details. 

2802.7 A doublet, places them at redshifts 2; > 1.246 and 2; > 1.238, re- 
spectively (Landt, 2012). PKS 0447-439 was detected as a very-high energy 
7— ray source with H.E.S.S., showing an energy spectrum from 200 GeV - 
1.5 TeV (Zech et al., 2011). 

The optical depth of 1 TeV photons at such large redshift is about 20 (see 
Figures 12 and 13), so that any flux of ~ 1 TeV photons is attenuated by a 
factor of ~ 10“®. This would suggest that the intrinsic 1 TeV luminosity of 
this blazar is about 10^° larger than that of Mrk 421! Any attempts to solve 
this problem by lowering the 7 — 7 opacity of the intervening intergalactic 
medium will require the unrealistic reduction of the EBL intensity below the 
lower limits determined by the IGL (see Figure 7). 

It is therefore unlikely that primary 7— rays could have reached the Earth 
from such distance, suggesting that most of the 7— rays from this object must 
be secondary photons created by the interaction of cosmic-ray protons at 
relatively close distances to the observer (Kusenko, 2012; Aharonian et ah, 
2012). The secondary photons are generated by the interactions of the pro- 
tons with the GMB and the EBL. For the secondary photons to be detected, 
they have to be produced within a distance A.y, the mean-free-path of the 
secondary 7— rays, of the Earth, and be deflected into the viewing angle of 
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the telescope (e.g Stanev et ah, 2000). For these photons to account for part 
of the energy spectra in blazars, their angular deviation from the primary 
source direction has to be smaller than the point spread function of the tele- 
scope (~ 0.1°). This requires the protons to travel on a ’’straight” trajectory 
until the interaction region, which sets a range 10 — 10 G on the 

line-of-sight intensity of the intergalactic magnetic held (Aharonian et ah, 
2012). While this value is signihcantly smaller than commonly accepted up- 
per limits for the intergalactic magnetic held (Kronberg, 1994, 2010), weak 
helds are not ruled out. 

Generally, leptonic models (SSG and EG) are capable of explaining the 
spectral energy distributions of most blazars. However, hadronic particle ac- 
celeration lies at the origin of the cascade model, potentially holds the key 
to understanding blazar jets, and may have far reaching consequences for 
their role in high energy astrophysics. Hadronic jet models are attractive for 
explaining the origin of ultra-high energy (UHE: E > 10^^ eV) cosmic rays 
that are detected by air shower arrays such as AUGER and HiRes (Kam- 
pert, 2012; Sokolsky & for the HiRes Gollaboration, 2010). Hadronic jet 
emission models were brought forward (Mannheim, 1993) around the time of 
the hrst detection of TeV photons from a blazar (Punch et ah, 1992). Pro- 
tons at energies of 10^® — 10^® eV are capable of reaching the threshold for 
photopion production with ambient photons, and generate a subsequent cas- 
cade inside the jet, and thereby make a substantial contribution to primary 
7— ray emission from blazars. Furthermore, proton synchrotron radiation 
and synchrotron emission from secondary particles contribute to the 7— ray 
spectrum. 

Those 7— ray emission models have been put to test through TeV 7— ray 
observations of nearby blazars for which EBL absorption is negligible. The 
acceleration of UHE protons along relativistic jets requires large magnetic 
helds, typically several 10s of Gauss. Upper limits to the size of the 7— ray 
emission region can be derived from the observation of TeV hares. The hrst 
big hare reported from Mrk 421 exhibited sub-hour scale hux variability (Gai- 
dos et ah, 1996). Assuming a Doppler factor of 10, the short hux variations 
set a limit to the size of the gamma-ray emission region (< 10^^ m) in the 
comoving frame of the relativistic jet (Gaidos et ah, 1996). Following Dermer 
(2012), a magnetic held strength of at least ~ 50 G is required to prevent 
10^® eV ions and protons from escaping the small emission region. Such large 
magnetic helds are still consistent with the substantial but plausible power 
requirements for the magnetic held energy (less than 0.1% of the Eddington 
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luminosity). Therefore, the observation of short flares does not pose a major 
obstacle to the viability of hadronic models. 

A consequence of hadron induced pair cascades in intergalactic space is a 
contribution to the diffuse 7 — ray background (DGB). Recent measurements 
of the DGB by Fermi (Abdo et al., 2010d) extent up to 100 GeV. Further 
measurements extending up to 1 TeV would provide important feedback to 
the presence of secondary 7 — rays from blazars. If indeed the intergalactic 
magnetic held strength is small, protons can travel cosmological distances 
and thus dissipate signihcant amounts of non-thermal energy into secondary 
7 — rays locally, and enhance the DGB well into the TeV regime. If the 
magnetic helds are relatively strong, distant sources would not contribute 
much to pair cascades initiated within the observer’s transparency zone, and 
one would expect the spectrum of the DGB to fall steeply above several 10s 
of GeV. 

However, the need to resort to hadronic acceleration model to explain 
the observed hux of PKS 0447-439 is probably premature. The redshift 
determination to this blazar has recently been challenged by Fumagalli et al. 
(2012) who independently obtained the PKS 0447-439 spectrum using the 
6.5 m Magellan Telescopes with a high S/N ratio of ~ 150 in the 6270- 
6300 A region of interest. They point out the existence of an atmospheric 
telluric absorption line at the wavelength of the claimed Mg II absorption 
line, which they showed was also present in the spectrum of two standard 
stars. As a result the redshift of PKS 0447-439 is still undetermined. 

9. EBL Limits on ’’Exotic” Energy Releases in the Universe 

The intensity and spectral shape of the EBL contain the memory of 
all energy releases in the universe since the epoch of recombination. The 
GORE/DIRBE limits on the UV to near-IR regions of the EBL proved use- 
ful in ruling out various ’’exotic” sources of energy in the early universe, such 
as decaying particles, exploding stars, or very massive objects (Dwek et al., 
1998, and references therein). 

Primordial (Population III) stars were suggested by Salvaterra & Fer- 
rara (2003) as the source of the excess 1-5 pm diffuse emission above the 
IGL intensity detected by Matsumoto et al. (2005). However, (Dwek et al., 
2005b) ruled out an extragalactic origin for the excess emission, since it would 
have produced a physically unrealistic intrinsic 7 — ray spectrum of the blazar 
PKS 2155-304. Furthermore, Dwek et al. (2005a) showed that such origin 


41 






Figure 15: The contribution of dark stars to the EBL for two different Limits on the EBL 
intensity. Lower limits, blue squares) are determined by the intensity of the IGL. Upper 
limits (red circles) are determined by absolute measurements of the EBL. The data used 
in the figure are listed in Tables 3-5 in bold. The shaded area depicts the range of the 
allowed EBL intensity as determines by UV to sub millimeter observations. 

would have required a Pop III star formation rate and energy output to be 
significantly higher than that predicted by hierarchical models for structure 
formation in a ACDM universe (Bromm & Loeb, 2002). Using theoretical 
limits on their formation rate, (Dwek et ah, 2005a) concluded that Pop III 
stars can contribute only a fraction of the EBL intensity. 

A more detailed study of the contribution of Pop III stars to the EBL 
intensity was conducted by Raue et al. (2009). Using 7 — ray derived con- 
straints on the EBL intensity they set a limit of 0.3 to 3 Mq Mpc“^ yr“^ on 
the formation rate of Pop III stars in the z = 7 — 14 redshift interval. 

More recently, (Maurer et ah, 2012) used EBL limits to constrain the 
properties of dark stars (DS) in the early universe. Dark stars are objects 
that have either accreted or captured weakly interacting massive particles 
(WIMPs, a dark matter candidate), which by annihilating inject energy into 
the stars before their radiative output is dominated by standard nuclear 
fusion processes. The formation rate of these stars, their luminosity and 
spectrum, and their effective lifetime are all free parameters of the model. 
Figure 15 depicts the EBL, to which the contribution of dark stars was added 
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to the IGL intensity, for two different sets of parameters. The blue curve rep- 
resents the EBL with a contribution of 106 Mq DS with a surface temperature 
of 5,000 K, and the red curve that with the added contribution of 690 M© DS 
with a surface temperature of 7,500 K. The colder dark stars are obviously 
ruled out, but the hotter ones are marginally consistent with current EBL 
limits, dominating the intensity of the IGL in the ~ 2 — 10 /rm wavelength 
region. 

10. Summary and Future Direction 

Very high energy 7— rays emitted from extragalactic sources are attenu- 
ated en route to earth by 7 — 7 interaction with EBL photons. 7— ray obser- 
vations can therefore be used to set limits on the EBL intensity, provided that 
the intrinsic 7— ray spectrum of the sources is known. Gonversely, knowledge 
of the EBL can be used to determine the intrinsic 7— ray spectrum of the 
different sources, thereby provide important constraints on mechanisms for 
their production. The main issues and results discussed in this review can 
be briefly summarized as follows: 

1. The EBL spectrum consists of two broad peaks, one at A ~ 1 /rm, 
representing the cumulative gravitational and nuclear energy releases 
by stars and AGNs over cosmic history. Their energy has been partially 
absorbed by dust and reradiated at IR wavelengths. This thermal 
dust emission component generates a second peak at ~ 100 — 200 /rm. 
Gurrent limits and detections of the EBL were presented in Figure 7, 
and the partitioning of its total intensity into the different emission 
components was presented in Table 6; 

2. The cross section for the 7—7 interaction is broad and peaks at energies 
E^iTeV) ~ 0.86 A(/r m). Gonsequently, the energy dependence of the 
7— ray opacity, a product of the cross section with the number density of 
EBL photons, reflects the spectral variation in the EBL. We identihed 
three potential breaks in the spectrum of 7— ray sources, one occurring 
between 10 and 500 GeV, a second break at 1 TeV and third at 10 TeV 
(see Figures 3 and 4); 

3. The spectral breaks at GeV and TeV energies have been used to set 
limits on the EBL intensity at near-IR wavelengths, and on the relative 
intensities of the peak of the stellar emission and the trough between 
the stellar and dust emission components of the EBL; 
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4. 7~ray derived limits on the EBL vary as different studies used different 
assumptions on the intrinsic 7— ray source spectra (see Figure 10). The 
strictest limits on the EBL are around the 1 and ~ 10 /rm wavelength 
regions of the EBL, with some approaching the lower limit on the EBL 
intensity set by the IGL (see Figure 11); 

5. Recent IR space and ground-based observations have resulted in closer 
agreement between the EBL limits and detections derived from mea- 
surements of the absolute sky brightness and lower limits set by the 
integrated light from galaxies. A summary of the recent observational 
status of the EBL is presented in Tables 3-5. The major gap in our 
knowledge of the EBL is in the A ~ 10—70 /rm wavelength region, where 
the thermal emission from interplanetary dust (the zodiacal light) dom- 
inates the brightness of the sky; 

6. Models of the EBL, employing various methods for determining the 
evolution of the galaxies’ spectral energy density with redshift seem 
to agree on the intensity and spectral shape of the UV to near-IR 
component of the EBL. However, there are considerable differences in 
their treatment of the redistribution of the intrinsic stellar and AGN 
output at IR wavelengths (see Figure 9); 

7. The discovery of TeV blazars with alleged redshifts ^ 1 has revived 
an alternative model for the creation of these high energy photons. In 
this model hadronic jets produce a cascade of ~ TeV 7— ray photons en 
route to earth, circumventing the attenuation problem from such high 
redshift sources (see §8). However, the redshift of these sources has 
been disputed, so the need to resort to such hadronic model is highly 
premature; 

8. Dark matter has been invoked to postulate the existence of a new class 
of primordial stars, powered by dark matter annihilation, instead of 
nuclear fusion. The energy release from these so-called dark stars could, 
in principle, lead to an observable signature in the EBL spectrum (see 
Figure 15); and hnally 

9. Observations suggest that the universe is essentially transparent to 
7— rays with energies < 2 TeV up to ~ 0.2, and energies of < 400 GeV 
up to 2: ~ 0.4 (see Figure 14). 

The future prospects of intensihed EBL studies with 7— rays are promis- 
ing, especially when considering the progress made over the last 5 years 
through the operation of Fermi and the current generation of atmospheric 
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Cherenkov telescopes. These instrument yielded combined energy spectra for 
at least 3 dozen extragalactic sources that cover up to 5 orders of magnitude 
in energy. These data have provided a hrst glimpse of the transition region 
from a transparent universe at 1 GeV, to TeV energies at which the universe 
gradually turns opaque with increasing redshift. 

Furthermore, Fermi detected a large set of blazars (Ackermann et ah, 
2011) including 310 FSRQs, 395 BL Lacertae objects and 156 candidate 
blazars, raising the specter for EBL studies with sizable source samples and 
different sources classes with CTA. Large samples of blazars have the po- 
tential to constrain the EBL in the optical/near-IR and mid- IR through a 
better understanding of the blazar subclasses (FSRQs, LBL, IBL, HBL) and 
their intrinsic spectra, and better photon statistics for the measurement of 
the redshift dependence of any spectral feature attributable to the EBL. 

Recent discoveries of TeV emission from nearby radio- and starburst 
galaxies also gives rise to future prospects of extending the reach of 7 — ray 
constraints up to the far-IR through energy spectra spanning a few GeV up 
to 100 TeV. These will be important to extend measurements of the 7 — ray 
opacity imposed by the EBL all the way from UV/optical/near-IR/mid-IR 
to the far-IR and thus provide additional constraints to the intensity ratios 
between the different wavelength regimes of the EBL. 

Limitations of precision EBL studies with 7 — rays arise from technical 
reasons, but are not unsurmountable. First, systematic uncertainties in the 
measurements of 7 — ray spectra with atmospheric Cherenkov telescopes arise 
from uncertainties in modeling the Earth’s atmosphere, which translates into 
uncertainties in the attenuation of Cherenkov light from electromagnetic cas- 
cades. To hrst order this affects the absolute energy scale of the measured 
7 — ray energies and is typically quoted at a level of 15% to 20%. Effects 
on the spectral index and shape are generally of second order but require 
detailed studies. Additional uncertainties lie in the absolute instrument cal- 
ibration of Cherenkov telescopes, e.g., mirror rehectivity, light losses in the 
focal plane, and uncertainties in the quantum efficiency of the photodetec- 
tors, affecting the light throughput and absolute energy scale. Uncertainties 
in the spectral indices and spectral shape are typically quoted at the level of 
AT = 0.05 — 0.1 or better, and depend mostly on the 7 — ray selection effi- 
ciency derived from Monte Carlo simulations and detailed detector modeling. 
Currently most energy spectra published (except for few haring sources) are 
dominated by statistical uncertainties and leave much room for improvement 
through a more sensitive instrument such as CTA with a collection area of 
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~ km?. At the same time calibration techniques are continually improving 
to provide a better handle on systematics as well. 

Astrophysical challenges facing EBL studies with 7— rays arise from the 
fact that currently only 55% of all Fermi detected BL Lacs have reliable red- 
shift measurements (Ackermann et ah, 2011). Attempts to increase the red- 
shift ident ideations through dedicated optical follow-up observations of the 
host galaxies of Fermi detected blazars during low flaring states are promis- 
ing to increase the fraction of blazars with known redshift. The redshift 
distribution for FSRQs detected by Fermi peaks at z = 1 extending to z = 
3.1, while BL Lacs peak at z = 0.2 reaching up to a redshift of z = 1.5. 
As already indicated, GeV/TeV detections of FSRQs are promising to allow 
constraints to the EBL in the UV/optical. Furthermore, the first detection 
of FSRQs by atmospheric Cherenkov telescopes out to redshifts oi z = 0.5 
provides the potential to complement EBL studies with a different class of 
sources, with different intrinsic spectral properties than BL Lacs and, most 
importantly, whose redshifts are known. 

Detection and limits on the EBL have been obtained by direct absolute 
measurements, and by galaxy number counts that provided lower limits on 
its intensity. Direct measurements are aided by ground-based or space-based 
observations that resolve the foreground emission from Galactic stars. Fu- 
ture observations can considerably improve direct absolute measurement of 
the EBL by determining the absolute brightness of the zodiacal light from 
high resolution observations of reflected solar Fraunhofer lines, or by making 
absolute sky measurements from the outer solar system, thus effectively re- 
moving any foreground contribution from the zodiacal cloud. Lower limits on 
the EBL can converge to the EBL itself with observations of sufficient depth 
and resolution to resolve all the galaxies that contribute to its intensity. 

The James Webb Space Telescope (JWST) is a large (6.6 m) IR space ob- 
servatory, passively cooled to temperatures below 50 K, that will be launched 
into orbit at the second Earth-Sun Lagrange point (L2). The wide band ob- 
servations with NIRGam at 2 - 5 /rm and MIRI at 6 - 25 /im will hll in crucial 
gaps in our knowledge of the EBL intensity at these wavelengths. NIRGam 
will cover a 10 aremin^ held of view (FOV) and will resolve galaxies down 
to a lOa hux limit of 11 nJy at 2 /xm in 10,000 s, making it about 10^ more 
sensitive than the Spitzer/IRAG 3.6 /rm hlter. With unparalleled resolu- 
tion and sensitivity, the JWST (Gardner et ah, 2006) will resolve the EBL 
at near-IR wavelengths. The MIRI instrument operating in its broadband 
imaging mode at 10 and 21 /xm will have a smaller FOV of ~2.6 aremin^. 
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and sensitivities of 700 nJy and 9 /iJy, respectively. It will resolve galaxies 
down to the confusion limit which, because of its large telescope size, will be 
signihcantly fainter than that of the Spitzer/MIPS instrument at 24 pm. 

Absolute measurements of solar absorption line prohles in the zodiacal 
light, for example the Fraunhofer Mg I line at 1.182819 pm or the Si I line 
at 1.210354 pm can be used to determine and remove the contribution of 
the ZL to the total sky brightness in the J (1.25 pm) band (Kutyrev et ah, 
2004, 2008). High resolution (R ~ 20,000) observations of these lines us- 
ing the ground-based spectrometer ZEFIR (Zodiacal Emission determination 
through Fraunhofer IR lines) will determine their equivalent width (EW) in 
the ZL spectrum. The continuum contribution of the ZL in the J-band can 
then be determined from knowledge of the EW of these lines in the solar 
spectrum. The DIRBE instrument provided absolute measurements of the 
sky brightness in this band, and ground-based sky surveys (2MASS) have re- 
moved most of the contribution of Galactic stellar emission from this band. 
Measurements of the absolute brightness of the ZL will thus provide the 
absolute intensity of the EBL in the DIRBE 1.25 pm waveband. Similar 
measurement could also be made from the ground in the 2.2 pm band. 

Finally, the contribution of the ZL to the observed sky brightness can be 
largely eliminated by mapping the absolute intensity of the diffuse emission 
from the sky from beyond the zodiacal cloud. A satellite mission conducting 
such EBL measurements from the outer solar system will be an important 
step in that direction (Cooray, 2011). 

The next decade will see signihcant improvement in y-ray technology, and 
measurements of the EBL, providing new insights into the spectrum of the 
different y-ray sources and the mechanism that generate them, and into the 
spectrum of the EBL which will provide new constraints on the history of 
all nuclear and gravitational energy releases in the universe. The synergy 
between TeV y-ray astronomy and EBL research will ensure that any devel- 
opment in each field will greatly benefit both. 
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Table 1: Glossary of Abbreviations of Spacecrafts^, Telescopes, and Instruments 
Abbreviation Full Name 

Akari Infrared imaging satellite (ASTRO-F) 

BLAST Balloon-borne Large-Aperture Submillimeter Telescope 

COBE Cosmic Background Explorer 

DIRBE Diffuse Infrared Background Experiment 

FIRAS Far Infrared Absolute Photometer 

CTIO Cerro Tololo Inter-American Observatory 

GALEX Galaxy Evolution Explorer 

Herschel Herschel Space Observatory 

PACS Photodetector Array Camera 

SPIRE Spectral and Photometric Imaging Receiver 

HST Hubble Space Telescope 

WFPC2 Wide Field Planetary Camera 

NICMOS Near IR Camera and Multi-Object Spectrometer 

IRTS Infrared Telescope in Space 

ISO Infrared Space Observatory 

ISOCAM ISO Camera 

JCMT James Clerk Maxwell Telescope 

SCUBA Submillimeter Common User Bolometer Array 

NTT New Technology Telescope 

Spitzer Spitzer Space Telescope 

IRAC Infrared Array Camera 

MIPS Multiband Imaging Photometer 

Subaru Optical, near-IR telescope 

Pioneer Interplanetary spacecraft 

WMAP Wilkinson Microwave Anisotropy Probe 

2MASS Two Micron All Sky Survey 

CTA Cherenkov Telescope Array 

Eermi Fermi Gamma-Ray Space Telescope 

H.E.S.S. High Energy Stereoscopic System 

lACT Imaging Air Cherenkov Telescope 

MAGIC Major Atmospheric Gamma-Ray Imaging Cherenkov Telescope 

Milagro Gamma-ray and cosmic-ray telescope 

VERITAS Very Energetic Radiation Imaging Telescope Array 

^ Spacecraft names are presented in italics, and their instruments are in- 
dented. 







Table 2: Extragalactic 7 — ray sources with GeV and TeV spectral information. 


Name 

Class 

redshift 

^GeV 

^TeV 

Range [TeV] 

References 

Ccntaurus A 

Radio 

0.0008 

2.76+0.05 

2. 7+0.5 

0.2 - 5 

[1], 12| 

M82 

SB 

0.00085 

2.2+0.2 

2. 5+0.6 

0.7 - 4 

[3], [4] 

NGC253 

SB 

0.00093 

1.95+0.4 

2.24+0.14 

0.3 - 50 

[3], [5] 

M87 

Radio 

0.0036 

2.17+0.07 

2.5+0.2 

0.2 - 10 

[6], [7], 8], [9] 

NGC 1275 

Radio 

0.018 

2.00+0.02 

3.96+0.37 

0.1 - 0.3 

[10], [11] 

IG 310 

Radio 

0.0188 

2.10+0.19 

2.0+0.14 

0.1 - 7 

[12], [13], [14] 

Markarian 421 

HBL 

0.031 

1.77+0.01 

2.48+0.03* 

0.1 - 5 

[15] 

Markarian 501 

HBL 

0.034 

1.74+0.03 

2.51+0.05^ 

0.1 - 10 

[16] 

lES 2344+514 

HBL 

0.044 

1.72+0.08 

2.78+0.09^ 

0.3 - 2 

[6], [17] 

Markarian 180 

HBL 

0.046 

1.74+0.08 

3.3+0.70 

0.2 - 1 

[6] [18] 

lES 1959+650 

HBL 

0.047 

1.94+0.03 

2.72+0.14 

0.2 - 2 

[6], [19] 

AP Lib* 

LBL 

0.048 

2.05+0.04 

2.5+0.2 

0.3 - 2 

[6], [20] 

BL Lacertae 

LBL 

0.069 

2.11+0.04 

3.6+0.5 

0.2 - 1 

[6], [21] 

PKS 2005-489 

HBL 

0.071 

1.78+0.05 

4.0+0.4 

0.2 - 2 

[6], [22] 

W Comae 

IBL 

0.103 

2.02+0.03 

3.81+0.35 

0.3 - 1 

[6], [23] 

PKS 2155-304 

HBL 

0.116 

1.84+0.02 

3.53+0.05 

0.4 - 5 

[6], [24] 

B3 2247+381 

HBL 

0.119 

1.84+0.11 

3.2+0.5 

0.2 - 1 

[6], [25] 

RGB J0710+591 

HBL 

0.125 

1.53+0.12 

2.69+0.26 

0.3 - 4.6 

[6], [26] 

H 1426+428 

HBL 

0.129 

1.32+0.12 

3.50+0.35 

0.3 - 10 

[6], [27] 

lES 1215+303 

IBL 

O.IS*^ 

2.02+0.02 

2.99+0.15 

0.1 - 1 

[6], [28] 

lES 0806+524 

HBL 

0.137 

1.94+0.06 

3.6+1. 0 

0.3 - 0.7 

[6], [29] 

IRXS J101015.9-311909 

HBL 

0.143 

2.24+0.14 

3.14+0.53 

0.3 - 1 

[6], [20] 

lES 1440+122 

IBL 

0.163 

1.41+0.18 

3.3+0.7 

0.3 - 1 

[6], [30] 

H 2356-309 

HBL 

0.165 

1.89+0.17 

3.09+0.24 

0.3 - 2 

[6], [31] 

VER J0648+152 

HBL 

0.179 

1.74+0.11 

4.4+0.8 

0.3 - 0.8 

[6], [32] 

lES 1218+304 

HBL 

0.184 

1.71+0.07 

3.07+0.09 

0.2 - 2 

[6], [33] 

lES 1101-232 

HBL 

0.186 

1.80+0.21 

2.88+0.17 

0.16 - 3.3 

[6], [31] 

RBS 0413 

HBL 

0.19 

1.55+0.11 

3.18+0.68 

0.25 - 1 

[6], [33] 

PKS-0447-439 

HBL 

0.205 

1.86+0.02 

4.36+0.49 

0.25 - 1 

[6], [34] 

lES 1011+496 

HBL 

0.212 

1.72+0.04 

4.0+0.50 

0.25 - 0.6 

[6], [35] 

lES 0414+009 

HBL 

0.287 

1.98+0.16 

3.44+0.27 

0.25 - 1.2 

[6], [36] 

S5 0716+714 

LBL 

0.31 

2.01+0.02 

3.45+0.54 

0.25 - 1.2 

[6], [37] 

lES 0502+675 

HBL 

0.416* 

1.49+0.07 

3.92+0.35 

0.25 - 1 

[6], [38] 

4G 21.35 

FSRQ 

0.43 

2.12+0.02 

3.75+0.27 

0.07 - 0.4 

[6], [39] 

3G 66A 

IBL 

0.44* 

1.85+0.02 

4.1+0.4 

0.22 - 0.45 

[6], [40] 

3G 279 

FSRQ 

0.536 

2.22+0.02 

3.03+0.9 

0.1 - 0.35 

[6], [41] 


★ Spectrum is not well fit by a powerlaw. A Spectrum shows variations. ^ redshift uncertain. ^ redshift 
recently given by Abdo et al. arXiV:1108.1420vl, different redshift was considered viable by Colin et al. 
2011 

References: [IJAbdo et al. (2010c); [2] Aharonian et al. (2009); [3] Abdo et al. (2010b); [4] Acciari et al. 
(2009b); [5] (Abramowski et al., 2012, for the H.E.S.S. collaboration); [6] Ackermann et al. (2011); [7] 
Berger et al. (2011a); [8] Acciari et al. (2010b); [9] Aharonian et al. (2006b); [10] Abdo et al. (2009); [11] 
Hildebrand et al. (2011); [12] Neronov et al. (2010) ; [13] Aleksic et al. (2010b); [14] Ackermann et al. 
(2011); [15] Abdo et al. (2011b); [16] Abdo et al. (2011a); [17] Acciari et al. (2011b); [18] Albert et al. 
(2006a); [19] Albert et al. (2006b); [20]Cerruti, M. for the H.E.S.S. Collaboration (2011); [21] Albert 
et al. (2007a); [22] Aharonian et al. (2005a); [23] Acciari et al. (2008); [24] Aharonian et al. (2007b); [25] 
Berger et al. (2011b); [26] Acciari et al. (2010a); [27] Petry et al. (2002); [28] Colin et al. (2011); [29] 
Acciari et al. (2009a); [30] Benbow (2011); [31] Aharonian et al. (2006a); [32] Errando, M. et al. (2011) 
[33] Acciari et al. (2010c); [33] §entiirk et al. (2011); [34] Zech et al. (2011); [35] Albert et al. (2007b); 
[36] Volpe et al. (2011); [37] Anderhub et al. (2009); [38] Benbow, W. et al. (2011); [39] Aleksic et al. 
(2011a); [40] Acciari et al. (2009d); [41] Aleksic et al. (2011b). 
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Table 3: Limits and detection of the extragalactic background light (EBL). 


A 


u m ^ sr 


Comment 

Reference 

(/xm) 

IGL 

EXT 

ABS 



0.1530 

0.68 ±0.10 

1.03 ± 0.15 


Galex 

[1] 

0.1595 

3.75 ± 1.25 



HST/STIS 

[2] 

0.2 

0.6 



FOCA/baiioon 

[3] 

0.2310 

0.99 ±0.15 

2.25 ± 0.32 


Galex 

[1] 

0.2365 




HST/WYVC2 

[2] 

0.30 



18± 12 

HST/WYVC2 

[4] 


2.7 ±0.3 

3.7 ± 0.7 


HST+gvonnA 

[5] 

0.36 

2.9t°;6 



HST+gvomiA 

[6] 

0.40 



< 36 (26 ± 10) 

dark cioud 

[7] 

0.44 

4.6t°o;5 


7.9±4.0 

Pioneer 10/11 

[8] 

0.45 



17S'T±ground 

[6] 


4.4 ±0.4 

6.1 ± 1.8 


17S'T±ground 

[5] 

0.5115 



< 39 (30 ±9) 

ground 

[9] 

0.55 



55 ±27 

HST/WYVC2 

[4] 

0.61 

6.0 ±0.6 

7.4 ± 1.5 


HST+gYO\xnd 

[5] 

0.64 



7.7±5.8 

Pioneer 10/11 

[8] 

0.67 




17S'T±ground 

[6] 

0.81 

8-oA^ 



17S'T±ground 

[6] 


8.1 ±0.8 

9.3 ± 1.6 


17S'T±ground 

[5] 

0.814 



57 ±32 

175'T/WFPC2 

[4] 

1.1 

Q 7+3.0 
^ -1.9 



17S'T±ground 

[6] 

1.25 



21±15 

COSE/DIRBE 

[10] 




54± 17 

CORE/DIRBE 

[11] 


10.9 ± 1.1 

11.5 ± 1.3 


17S'T±ground 

[5] 


1 1 7+5-6 



Subaru 

[12] 

1.4-4 



~ 60 - 15 

IRTS 

[13] 

1.6 

9.0+11 



HST+gvo\mA 

[6] 


n.5+_n 



Subaru 

[12] 

2.12 

10.0+1 



Subaru 

[12] 

2.2 

7.9t+ 



17S'T±ground 

[6] 


8.3 ±0.8 

9.0 ± 1.2 


17S'T±ground 

[5] 




20±6 

CORE/DIRBE 

[10] 




28 ±7 

CORE/DIRBE 

[11] 

3.5 



13.3 ±2.8 

CORE/DIRBE 

[10] 




13.8±3 

CORE/DIRBE 

[15] 

3.6 

5.4 

Q 0^^ ^ 


Spitzer/lYlAC 

[16] 




Spitzer/lYlAG 

[14] 

4.5 

3.5 



Spitzer/lYlAC 

[16] 

4.9 



22±12 

CORE/DIRBE 

[17] 

5.8 

3.6 



Spitzer/IRAC 

[16] 

8.0 

2.6 



Spitzer/lYlAC 

[16] 


ABS=absolute measurement; IGL=integrated galactic light; STK=lower limits from stacking analysis; 
EXT=extrapolated intensity from dN/dS 

^ Calculated for a 2.2 /rm intensity of 20.0 nW m~^ sr“^ and the Kelsall et al. (1998) ZL model. 


References: [1] Xu et al. (2005); [2] Gardner et al. (2000); [3] Milliard et al. (1992); [4] Bernstein (2007); [5] 
Totani et al. (2001); [6] Madau Sz Pozzetti (2000); [7] Mattila (1990),Leinert et al. (1998); [8] Matsuoka 
et al. (2011); [9] Dube et al. (1979),Leinert et al. (1998); [10] Levenson et al. (2007); [11] Cambresy et al. 
(2001); [12] Keenan et al. (2010); [13] Matsumoto (2001); [14] Levenson & Wright (2008); [15] Dwek Sz 
Arendt (1998); [16] Fazio et al. (2004); [17] Arendt Sz Dwek (2003) 


51 





Table 4: Limits and detection of the extragalactic background light (EBL). 




vh (nW 

m ^ sr 




A(^m) 

IGL 

STK 

EXT 

ABS 

Comment 

Reference 

15 

2.4 ±0.5 
2 . 7 ± 0.6 

1.9 ±0.5 




reO/ISOCAM 

KO/ISOCAM 

Akari 

Elbaz et al. (2002) 
Metcalfe et al. (2003) 
Hopwood et al. (2010) 

16 

2.2 ±0.2 




Spitzer 

Teplitz et al. (2011) 

24 

1.9 ±0.6 
1.8 ±0.2 
2.29 ± 0.09 


^•'- 0.7 
2.0 ± 0.2 

2 . 86 t °;« 


Spitzer/MlPS 

Spitzer/MlPS 

Spitzer/MlPS 

Papovich et al. (2004) 
Chary et al. (2004) 
Bethemiin et al. (2010a) 

60 




28.1±1.S±7 

COBE/DIKBE 

Finkbeiner et al. (2000) 

65 




12.5±1.4±9.2 

Akari 

Matsuura et al. (2011) 

70 

5.4 ±0.4 
4.52 ± 1.18 

7.1 ± 1.0 

7.4 ± 1.9 
6 6^’^''^ 


Spitzer/MlPS 

Herschel/PACS 

Spitzer/MlPS 

Herschel/PACS 

Prayer et al. (2006) 
Dole et al. (2006) 
Bethermin et al. (2010a) 
Berta et al. (2011) 

90 




22.3±1.7±4.7 

Akari 

Matsuura et al. (2011) 

100 

8.35 ±0.95 

9 . 4 ± 1.1 


< 34 (22 ±6) 
12.5 ± 5 

24.6 ±2.5 ±8 
23.4 ±6.3 
14 . 4 ± 6.0 

6.6 ± 1.8 ±2.1 

COBE/DIRBE (D, KZL) 
COBE/EIKBE (D, WZL) 
COBE/DIRBE 
COBE/DIRBE (D, KZL) 
DIRBE (F, WZL) 
Herschel/PACS 
Spitzer/MlPS 

Hauser et al. (1998) 
Wright (2004) 
Finkbeiner et al. (2000) 
Lagache et al. (2000) 
Dole et al. (2006) 
Berta et al. (2010) 
Benin et al. (2012) 

140 




25.0 ±6.9 

15.0 ±5.9 
32± 13 

24.2 ± 11.6 
22 ±7 

20.1 ±3.4 ±1.1 

COBE/DIRBE (D, KZL) 
COBE/DIRBE (F, KZL) 
COBE/DIRBE 
COBE/DIRBE (D, KZL) 
COBE/DIRBE (D, WZL) 
Akari 

Hauser et al. (1998) 
Odegard et al. (2007) 
Schlegel et al. (1998) 
Lagache et al. (2000) 
Wright (2004) 
Matsuura et al. (2011) 





12.4 ±6.9 
12.6±6.0 

COBE/DIRBE (F, WZL) 
COBE/FIRAS 

Dole et al. (2006) 
Fixsen et al. (1998) 
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Table 5: Limits and detection of the extragalactic background light (EBL). 




(nW 

m ^ sr 




A(^m) 

IGL 

STK 

EXT 

ABS 

Comment 

Reference 

160 

9.49 ±0.59 
8.9± 1.1 

13.4 ± 1.7 

11.4 ±0.7 


13.7 ±3.9 ±0.8 

14.4 ±0.8 ±2.3 
13.7±6.1 

Herschel/'PAC'^ 

Akari 

Herschel/FACS 

Spitzer/'MWS 

Spitzer/yUPS 

COBE/FIKAS 

Dole et al. (2006) 
Matsuura et al. (2011) 
Berta et al. (2010) 
Bethermin et al. (2010a) 
Benin et al. (2012) 
Fixsen et al. (1998) 

170 

19.1 ±5.6 ±5.3 




ISOPHOT 

Juvela et al. (2009) 

240 




13.6 ±2.5 
13 ± 2.5 

12.7 ± 1.6 
17 ±4 

11.0 ±6.9 
12.3 ±2.5 
10.&±4.3 

COBE/DIRBE (D, KZL) 
COBE/DIRBE (D, WZL) 
COBB/DIRBE (F, KZL) 
COBB/DIRBE 
COBB/DIRBE (D, KZL) 
COBB/DIRBE (F, WZL) 
COBB/FIRAS 

Hauser et al. (1998) 
Wright (2004) 
Odegard et al. (2007) 
Schlegel et al. (1998) 
Lagache et al. (2000) 
Dole et al. (2006) 
Fixsen et al. (1998) 

250 

1.73 zb 0.33 
1.55 zb 0.30 

8.6 zb 0.6 
O.U_2 6 

7.40 zb 1.42 

10.13l|®° 

10.3±4.0 

BLAST 

BLAST 

Herschel/SPIKE 

Herschel/SPIKE 

FIRAS 

Marsden et al. (2009) 
Bethermin et al. (2010b) 
Oliver et al. (2010) 
Bethermin et al. (2012) 
Fixsen et al. (1998) 

350 

U.UD_q Q4 

0.63 zb 0.18 
0.77 zb 0.16 

4.93 ± 0.34 

9 

^•^-2.0 
4.50 zb 0.90 

6A6t\fr 

5.6±2.1 

BLAST 

BLAST 

Herschel/SPIKE 

Herschel/SPIKE 

FIRAS 

Marsden et al. (2009) 
Bethermin et al. (2010b) 
Oliver et al. (2010) 
Bethermin et al. (2012) 
Fixsen et al. (1998) 

500 

U.Ui_0.01 
0.15 zb 0.07 
0.14 zb 0.03 

2.27 zb 0.20 
1 4+2.1 

1.54 zb 0.34 

2.80t°j? 

2.4±0.9 

BLAST 

BLAST 

Herschel/SPIKE. 

Herschel/SPIKE 

FIRAS 

Marsden et al. (2009) 
Bethermin et al. (2010b) 
Oliver et al. (2010) 
Bethermin et al. (2012) 
Fixsen et al. (1998) 

850 

0.12 ±0.03 
0.24 ±0.03 



0.5±0.21 

SCUBA 

SCUBA 

COBB/FIRAS 

Coppin et al. (2006) 
Zemcov et al. (2010) 
Fixsen et al. (1998) 

200-1000 





FIRAS 

Fixsen et al. (1998)^ 

1 a = (1.3± 

= 0.4) X 10-5; 

k = 0.64 ±0.12; T = 

(18.5 ± 1.2) K-, Ao 

= 100 
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Table 6: Total EBL Intensity of Different Models (nW m ^ sr 


Model 

Stars^ 

Dust^ 

Total 

Franceschini et ah (2008) 

25 

40 

65 

Dominguez et al. (2011) 

25 

44 

69 

Gilmore et al. (2011) 

25 

23 

48 

Mazin & Raue (2007) 

30 

26 

56 

Finke et al. (2010) 

27 

20 

47 

Stecker et al. (2006) 

61 

35 

96 

Observations^ 

23-93 

20-110 

42-202 

From CSFR4 

- 

- 

21-66 


^ Integrated intensity from 0.1 to 10 /tm. 

^ Integrated in tensity from 10 to 1000 /rm. 

3 From limits and detections of the EBL (Figure 7) 

4 Total integrated intensity inferred from the CSFR (Fig. 8) 
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